Fundamental Theoretical Bias in Gravitational Wave Astrophysics 
and the Parameterized Post-Einsteinian Framework 
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We consider the concept of fundamental bias in gravitational wave astrophysics as the assump- 
tion that general relativity is the correct theory of gravity during the entire wave-generation and 
propagation regime. Such an assumption is valid in the weak field, as verified by precision experi- 
ments and observations, but it need not hold in the dynamical strong-field regime where tests are 
lacking. Fundamental bias can cause systematic errors in the detection and parameter estimation of 
signals, which can lead to a mischaracterization of the universe through incorrect inferences about 
source event rates and populations. We propose a remedy through the introduction of the pa- 
rameterized post-Einsteinian framework, which consists of the enhancement of waveform templates 
via the inclusion of post-Einsteinian parameters. These parameters would ostensibly be designed 
to interpolate between templates constructed in general relativity and well-motivated alternative 
theories of gravity, and also include extrapolations that follow sound theoretical principles, such as 
consistency with conservation laws and symmetries. As an example, we construct parameterized 
post-Einsteinian templates for the binary coalescence of equal-mass, non-spinning compact objects 
in a quasi-circular inspiral. The parametrized post-Einsteinian framework should allow matched 
filtered data to select a specific set of post-Einsteinian parameters without a priori assuming the 
validity of the former, thus either verifying general relativity or pointing to possible dynamical 
strong-field deviations. 

PACS numbers: 04.80.Cc,04.80.Nn,04.30.-w,04.50.Kd 



I. INTRODUCTION 

A. The concept of bias and 
the validity of General Relativity 

The problem of bias has had a negative impact on 
many scientific endeavors by causing scientists to draw 
incorrect inferences from valid data. As the next gen- 
eration of gravitational wave (GW) detectors [E S [| 
promises to open a new observational window onto the 
universe, one wonders what bias might be present in this 
new branch of astronomy that could lead to incorrect 
conclusions about the universe seen in GWs. This is a 
particularly pertinent question to address early on, as the 
first observations are expected to have quite low signal- 
to-noise ratio (SNR). 

Bias is here to be understood as some set of a pri- 
ori assumptions, prejudices or preconceptions that affect 
conclusions derived from properly collected data. We 
purposely exclude from this definition experimentally- 
induced systematic errors, as we wish to concentrate on 
errors introduced during post-processing and analysis. In 
GW astrophysics, a principal source of bias is the a pri- 
ori assumption, often unstated, that general relativity 
(GR) is the correct theory that describes all gravitational 
phenomena at the scales of relevance to GW generation 
and propagation. This bias is ingrained in many of the 
detection and parameter estimation tools developed to 
mine GW data, and extract astrophysical information 
from them (see for example Q and references therein). 

The study of theoretical bias in GW astrophysics is not 
completely new. Recently, Cutler and Vallisneri Q ana- 



lyzed the issue of systematic error generated by the use of 
inaccurate template families. This issue can be broadly 
thought of as a modeling bias, where the preconception 
relates to physical assumptions to simplify the solutions 
considered (e.g. that all binaries have circularized prior 
to merger), or unverified assumptions about the accuracy 
of the solution used to model the given event. The study 
in Q considered inaccuracies arising from errors in solu- 
tions to the Einstein equations due to the use of the post- 
Newtonian (PN) approximation scheme, and they found 
these inaccuracies could dominate the error budget. This 
paper differs from that analysis in that we are concerned 
with the validity of the Einstein equations themselves, an 
assumption that can be considered a fundamental bias. 

Systematic errors created by fundamental bias may be 
as large, if not larger, than those induced by modeling 
bias, as waveforms could deviate from the GR prediction 
dramatically in the dynamical strong-field, if GR does 
not adequately describe the system in that region. This 
is particularly worrisome for template-based searches, as 
the event that will be ascribed to a detection will be 
the member of the template bank with the largest SNR. 
Given that GR is quite well tested in certain regimes, 
many sources cannot have deviations so far from GR 
as to prevent detection with GR templates (albeit with 
lower SNR). Thus, if templates are used based solely 
on GR models, although the corresponding events may 
be "heard" , any unexpected information the signals may 
contain about the nature of gravity will be filtered out. 

But is it sensible to expect GR deviations at all? GR 
has been tested to high accuracy in the Solar System and 
with binary pulsars [aQ- However, in the Solar System 
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gravitational fields are weak and particle velocities are 
small relative to the speed of light, and all such tests 
only probe linear perturbations in the metric beyond a 
flat Minkowski background. Binary pulsar tests begin to 
probe the strong-field regime in terms of the compact- 
ness of the source {i.e. the ratio of the source mass to 
its radius), though known systems are still quite weak 
with regard to the strength of the dynamical gravita- 
tional fields (characterized by the ratio of the total mass 
of the system to the binary separation). In this sense 
then, the dynamical, strong-field region of GR has so far 
eluded direct observational tests [1| . 

The lack of experimental verification of GR in the dy- 
namical strong-field could be remedied through GW ob- 
servations. For example, a compact object such as a 
black hole (BH) or neutron star (NS) spiraling into a 
supermassive BH will emit waves that carry a map of 
the gravitational field of the central BH, a program of 
relevance to the planned Laser Interferometer Space An- 
tenna mission (LI SA) M , and sometimes referred to as 
bothrodesy [ill . [l2 . [l3j . Such maps would then allow 
for detailed tests of alternative theories of gravity. Ad- 
ditional work on constraining alternative theories have 
concentrated on using either the binary inspiral or the 
post-merger ringdown phase. With the former, studies 
have focused on the search for the gray iton Compton 
wavelength [H, [H, [H, [H [H [H, HO, |21| . the existence 
of a scalar component to th e g ravitational interaction 
[13, [H, [H, iO, im, m m him and the existence 
of gravitational parity violation [26l [27j. Studies of the 
ringdown phase have concentrated on violations of the 
GR no-hair theorem (that the mass and spin completely 
determine the gravitational field of a rotating BH) by 
proposing to verify certain consistency relations between 
quasinormal ringdown (QNR) modes [H, [2^, [sO, HH . All 
of the above tests, however, have focused on specific al- 
ternative theories and have not investigated the issue of 
fundamental bias that we wish to address here. 



B. Toward a framework for detecting fundamental 
bias in GW observations of compact object mergers 

In this paper, we are not so much interested in specific 
tests of particwZar alternative theories, but rather we want 
to introduce a framework that may allow one to quantifi- 
ably investigate the consequences of fundamental bias in 
GW astronomy. We propose to do so via the introduc- 
tion of the parameterized post-Einsteinian (ppE) frame- 
work, an alog ous to the parameterized post-Newtonian 
(ppN) [U m, m, Isp . [3^ or the parameterized post- 
Keplerian (ppK) ones [a.[37l[38l|. where model waveforms 
are enhanced in a systematic and well-motivated manner 
by parameters that can measure deviations from GR. The 
construction of such a ppE framework is a quixotic task 
in general, and so to make it more manageable we will 
concentrate on a subset of possible GW sources: the in- 
spiral, merger and ringdown of binary JiH-like compact 



objects. 

Given that we want to liberate GW astronomy from 
the assumption that GR is correct in all regimes, what 
sense then does it make to begin with a binary BH 
merger, an event wholly within the realm of GR? First, 
there is strong observational evidence that highly com- 
pact objects exist, and in many cases the observations are 
consistent with the supposition that they are BHs [39| : 
the high luminosity of quasars and other active galactic 
nuclei (AGN) can be explained by gravitational bind- 
ing energy released through gas accretion onto super- 
massive (10^ — IO^Mq) BHs at the centers of the galax- 
ies [13, Ell ; several dozen X-ray binary systems discov- 
ered to date have compact members too massive to be 
NSs, and exhibit phenomena consistent with matter in- 
teractions originating in the strong-field regime of an in- 
ner accretion disk [44I; the dynamical motions of stars 
and gas about the centers of nearby galaxies and our 
Milky Way Galaxy suggest the presence of very mas- 
sive, compact objects, the most plausible explanation 
being supermassive BHs [i^ E3, E3- Therefore, it is 
not too bold to assume that even if GR does not accu- 
rately describe compact objects, whichever theory does 
must nevertheless still permit BH-Zifce solutions. Second, 
from studies of the orbital decay of the Hulse- Taylor (46j 
and more recently discovered binary pulsars [47l . E3 > 
is again a rather conservative conclusion that the early 
evolution of binary compact objects in the universe is ad- 
equately governed by GR — namely, that binary systems 
are unstable to the emission of quadrupole gravitational 
radiation. 

Based on the above considerations, starting with GR 
binary BH merger waveforms seems sound. We are then 
faced with the question of how to modify these waveform 
in a sensible manner. In theory, there are uncountably 
many conceivable modifications to GR that only manifest 
in the late stages of the merger. To make this question 
manageable, we shall guide our search for ppE expansions 
by looking to alternative theories that satisfy as many of 
the following criteria as possible: 

(i) Metric theories of gravity: theories where gravity 
is a manifestation of curved spacetime, described 
via a metric tensor, and which satisfies the weak- 
equivalence principle [ioj . 

(ii) Weak-field consistency: theories that reduce to GR 
sufficiently when gravitational fields are weak and 
velocities are small, i.e. to pass all precision, exper- 
imental and observational tests. 

(iii) Strong-field inconsistency: theories that modify GR 
in the dynamical strong-field by a sufficient amount 
to observably affect binary merger waveforms. 

Notice that the weak-field consistency criterion also re- 
quires the existence and stability of physical solutions, 
such as the Newtonian limit of the Schwarzschild met- 
ric to describe physics in the Solar System. One might 
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also wish that other criteria be satisfied, such as well- 
posedness of the initial-value problem, the existence of 
a well-defined, relativistic action, and that the theory 
be well-motivated from fundamental physics (eg. string 
theory, loop quantum gravity, the Standard Model of el- 
ementary interactions, etc), but we shall not impose such 
additional requirements here. 

Instead of concentrating on a particular theory that 
satisfies the above criteria, here we are more concerned 
with being able to measure generic deviations from GR 
predictions. In spite of the possibly uncountable num- 
ber of modifications one can introduce to the GR action, 
parameterizations of generic deviations is easier if one 
concentrates on the waveform observable, the response 
function, which in turn can be described by a complex 
function in the frequency domain. Deviations from GR 
then translate into deviations in the waveform amplitude 
and phase. Such deviations could be induced, for exam- 
ple, by modifications to the GW emission formulae (e.g. 
the quadrupole formula), new polarization modes, new 
propagating degrees of freedom and new effective forces. 
In order to ease the flow of the paper, we will not further 
discuss details of how certain alternative theories modify 
the gravitational waveform here, referring the interested 
reader to Appendix A. 



C. Two key questions related to fundamental bias 

In view of the above considerations, we shall attempt 
to lay the foundations to answer the following two ques- 
tions related to fundamental bias in GW astronomy, here 
focusing on binary compact object mergers: 

(i) Suppose gravity is described by a theory differing 
from GR in the dynamical, strong-field regime, but 
one observes a population of merger events filtered 
through a GR template bank. What kinds of sys- 
tematic errors and incorrect conclusions might be 
drawn about the nature of the compact object pop- 
ulation due to this fundamental bias? 

(ii) Given a set of observations of merger events ob- 
tained with a GR template bank, can one quantify 
or constrain the level of consistency of these obser- 
vations with GR as the underlying theory describing 
these events? 

As we will discuss in the next subsection, the ppE 
framework will be able to answer question (2), at least 
for the class of deviations from GR considered here. As 
for question (1), we shall not specifically address it in this 
paper, though will here give a couple of brief examples 
to clarify the question. As a first example, suppose the 
"true" theory of gravity differs from GR in that scalar 
radiation is produced in the very late stages of a merger, 
and during the ringdown. This will cause the inspiral 
to happen more quickly compared to GR, resulting in 
a late-time dephasing of the waveform relative to a GR 



template. Also, less power may be radiated in GWs dur- 
ing the ringdown phase. In all then, these events may be 
detected by GR templates, though with systematically 
lower SNR, i.e. be seemingly more distant. Here then, 
the fundamental bias could make one incorrectly infer 
that merger events were more frequent in the past. 

For a second hypothetical example, consider an ex- 
treme mass ratio merger, where a small compact object 
spirals into a supermassive BH. Suppose that a Chern- 
Simons (CS)-like correction is present, altering the near- 
horizon geometry of the BH as described in [s^, HH- 
To leading order, the CS correction reduces the effective 
gravitomagnetic force exerted by the BH on the compact 
object; in other words, the GW emission would be sim- 
ilar to a compact object spiraling into a GR Kerr BH, 
but with smaller spin parameter a. Suppose further that 
near-extremal (a w 1) BHs are common (how rapidly as- 
trophysical BHs can spin is an interesting and open ques- 
tion). Observation of a population of CS-modified Kerr 
BHs using GR templates would systematically underes- 
timate the BH spin, leading to the erroneous conclusion 
that near-extremal BHs are uncommon, which could fur- 
ther lead to incorrect inferences about astrophysical BH 
formation and growth mechanisms. 



D. Toward a ppE Construction 

The concept of a ppE framework is in close analogy 
to the ppN one, proposed by Nordtvedt and Will [1, [H, 
m. 111, H^, to deal with an outbreak of alternative 
theories of gravity in the 1970's. Such a framework al- 
lows for model-independent tests of GR in the Solar Sys- 
tem, through the introduction of ppN parameters in the 
weak-field expansion of the metric tensor. When these 
parameters take on a specific set of values, the metric 
tensor becomes identical to that predicted by GR, while 
when it takes on different values, the gravitational field 
becomes that predicted by certain alternative theories. 
In the same way, we introduce ppE parameters to inter- 
polate between different theories, but we parameterize 
the GW response function instead of the metric tensor, 
as the former is the observable in GW astrophysics. 

We shall follow the ppN route in the construction of a 
ppE framework: we will explore a set of alternative the- 
ories and their effects on GWs, and from these, we will 
phenomenologically infer and engineer a ppE template 
family that not only captures the waveforms from known 
theories, but also other phenomenological corrections. As 
we shall see, our approach reproduces corrections to the 
GW response functions predicted by a large class of al- 
ternative theories of gravity proposed to date, including 
Brans-Dicke theory, massive graviton theories and non- 
dynamical CS modified gravity. 

As a first step, we shall not consider here the effect of 
scalar, vectorial or all six tensorial polarizations on the 
response function, nor will we attempt to classify alter- 
native theories based on the different perturbative modes 
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that can be excited. In fact, such a classification already 
exists, the E{2) scheme 0, HI], which organizes differ- 
ent alternative theories according to whether they excite 
certain contractions of the Weyl tensor as characterized 
by the Newman-Pcnrose (NP) scalars. Here we will con- 
centrate on the two "plus" and "cross" tensorial degrees 
of freedom, and how these change in alternative theories. 
Additional excited modes will only be considered insofar 
as they affect the source dynamics, and thus indirectly 
the structure of the plus and cross modes. 

In the future it would be useful to extend the ppE 
templates to include additional GW polarizations and 
modes. However that would be a non-trivial effort, the 
main reason being that currently there is a lack of al- 
ternative theories that simultaneously satisfy the criteria 
discussed above and allow for the excitation of additional 
modes. For example, Brans-Dicke theory does possess 
a breathing mode, though it is strongly suppressed by 
constraints on its coupling constant. Similarly, Einstein- 
Aether theory generically has 5 propagating degrees of 
freedom, but these decouple and radiation can be shown 
to remain quadrupolar [53|. Thus, unlike with the other 
ppE extensions considered here, there is little theoretical 
guidance available on how GR waveforms would be mod- 
ified if additional GW polarizations are excited. Without 
such guidance then, it will be almost impossible to argue 
that any proposed extensions are well-motivated. 

A separate issue is whether the effect of additional 
polarization states can be observed with contemporary 
or planned detectors. For direct measurement one re- 
quires 2N arms to be sensitive to N polarization states 
[54| . for example using multiple ground-based detectors, 
Doppler tracking of multiple spacecraft, or pulsar timing 
arrays [11] . Alternatively, a space-borne detector such as 
LISA can be sensitive to multiple polarization states for 
events that last a sizable fraction of the spacecraft's or- 
bit, as the motion of detector then effectively samples the 
waveform with multiple arm orientations. For ground- 
based detectors, binary merger events happen too quickly 
for a similar strategy to be effective. However, other 
sources, for example the waves that could be emitted by 
a "mountain" on a pulsar, will produce coherent GWs 
over many rotation periods of the Earth, and thus could 
also contain information on multiple polarization states. 
We will not consider these interesting issues here. 



E. A ppE template family for the inspiral, merger 
and ringdown of black hole-like compact objects 

For this initial study, we construct ppE templates 
describing only the quasi- circular coalescence of non- 
spinning and equal-mass compact objects. Even within 
this restricted class of events, the ppE construction is 
non-unique, and certainly more refined versions could be 
developed. We introduce several ppE template families 
that vary in the number of ppE parameters, and thus, in 
the amount of fundamental bias each family assumes in 



its construction. The frequency-domain (overhead tilde) 
ppE templates with the least number of ppE parameters 
that we derive is the following (more general representa- 
tions are provided in Sec. lVCp : 

p|^"'(/).(l + au«)e*'3"' /</:m, 

h{f) = I 7U-e^(*+-) /:m < / < /mhd, (1) 

where the subscript IM and MRD stand for inspiral- 
merger and merger-ringdown. In the inspiral phase 
(/ < /im), the GW is described by a "chirping" complex 
exponential, consisting of the GR component h['^^\f) 
corrected by ppE amplitude and phase functions with 
parameters (a, a, (3, b). Here u = nMf is the inspiral re- 
duced frequency, M = Mrf"^^ is the chirp mass with 
symmetric mass ratio 77 = mim2/M'^ and total mass 
M = mi -|- 7712 (though again, we will only focus on the 
equal mass case). The merger phase (/im < / < /mrd) is 
treated as an interpolating region between inspiral and 
ringdown, where the merger parameters (7,(5) are set by 
continuity, and the merger ppE parameters are (c, e) . In 
the ringdown phase (/ > /mrd), the GW is described by a 
single-mode generalized Lorentzian, with real and imagi- 
nary dominant frequencies /^d and r, ringdown parame- 
ter ^ also set by continuity, and ppE ringdown parameters 
(k, d). The transition frequencies (/im,/mbd) can either 
be treated as ppE parameters, or (for example) set to 
the GR light-ring frequency and the fundamental ring- 
down frequency, respectively. In a later section, we shall 
present more general ppE waveforms with inherently less 
fundamental bias, some of which are reminiscent to the 
work of [H, m, [S^l , except that here we are interested 
in tests of GR through the determination of ppE param- 
eters, instead of the measurement of FN ones. 

The ppE template presented above is non-unique; in 
fact, the ppE framework as a whole is inherently non- 
unique, as are the ppN or ppK ones, because a finite 
parameterization cannot represent an infinite space of al- 
ternative theory templates. However, this ppE family is 
minimal within the class of templates considered here, 
as it employs the smallest number of ppE parameters 
necessary to reproduce corrections to the GW response 
function from well-known alternative theories of gravity 
in the inspiral phase: 

(i) GR is reproduced with {a,a,f3,b) — (0, a,0, 6), 
(c,e) = (-2/3,1) and («:, d) = (1, 2); 

(ii) Jordan-Brans-Dicke-Fierz, or simply Brans-Dicke 
theory (BD), with {a,a,P,b) = (0, a, /3bd, — 7/3), 
and Pbtd related to the Brans-Dicke coupling param- 
eter (see Eq. (KE\i ): 

(iii) Massive graviton (MG) theories with (a, a, (3, b) — 
(0, a, /3mg, — 1), and /3mg related to the graviton 
Compton wavelength (see Eq. (|A10p ): 
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(iv) CS modified gravity with (a, a, /?, b) = (ofcs, 1, 0, 6), 
and Qfcs related to the CS couphng parameter (see 
eg. Eq. (jA18|). 

The allowed range of the ppE parameters is not com- 
pletely free, as Solar System and binary pulsar experi- 
ments have already constrained some of them. The sub- 
set of inspiral ppE templates that reproduce well-known 
theories after fixing a and b is then constrained to some 
range in (a, (3). More precisely, since the Brans-Dicke 
coupling parameter has been constrained to lu-bo > 4x 10* 
by the Cassini spacecraft [s^, this automatically forces 
/3bd < 3.5xlO"**(si-S2)^?7^/^ when 6 = -7/3, where si,2 
are the sensitivities of the binary components. For ex- 
ample, for a IAMq binary NS system, P^d ^ 10~^ when 
b = —7/3. Similarly, since the Compton wavelength 
of a massive graviton has already been constrained to 
\g > 3.4 X 10^^ km from pulsar timing observations [HH, 
this implies /3mg < 8.5 x lO^^^DMil + z)-^ km"^ when 
6 = — 1, where D is a distance measure to the source. 
For example, for a 1.4M0 binary NS system at a red- 
shift of z = 0.1, one finds the prior /?mg ^ 1-6 x 10~^ 
when 6 = — 1. Finally, the non-dynamical CS coupling 
parameter has been constrained by binary pulsar obser- 
vations to be such that for an aligned binary (zero in- 
clination angle), acs ^ 70 x 10~'^'^£'/A^ when a = 1, 
which for a binary NS at redshift z = 0.1 translates to 
acs < 4 X 10-1". 

One can also consider placing priors on other ppE pa- 
rameters introduced in Eq. ([1]) without fixing the wave- 
form to represent a particular well-known theory. How- 
ever, since the merger and the ringdown have not been 
constrained at all by observation, one cannot really re- 
strict (c, e, K, c?). On physical grounds, one can only re- 
quire that d > and real and (e, k) also real, such that 
the waveform is well behaved at large frequencies. Since 
the quadrupole formula has been verified to leading or- 
der with binary pulsar observations, one also expects that 
for low- frequencies au" < 1, which for example implies 
that if < a ^ 1 then \a\ < 1, while if a 1 then 
a is essentially unconstrained. If one docs not expect 
gravitational radiation to be sourced below dipolar or- 
der one must have b > —7/3; similarly, if one expects 
strong-field, source generation deviations to only arise 
beyond quadrupolar order then b > —5/3. Of course, 
the stronger the priors imposed, the stronger the bias 
the ppE model inherits. 

F. The Possible Role of the ppE Framework in 
Data Analysis 

We do not propose here to employ the ppE templates 
for direct detection, but rather for post- detection analy- 
sis. Following the detection of a GW by a pure GR tem- 
plate, that segment of data could then be reanalyzed with 
the ppE templates. In principle, one should search over 
all system and ppE parameters, but in practice, since 
the ppE waveforms are deformations of GR templates. 



one can restrict attention to a neighborhood of the sys- 
tem parameter space centered around the best-fit ones 
obtained through filtering with a GR waveform. 

The templates presented above do contain some mod- 
eling error, i.e., the subset of the ppE template family 
with all parameters set to the GR values — the ppE-GR 
templates — are not the "exact" solution to the event in 
GR (for which there is no known closed form solution). 
What is more important though is that the ppE-GR tem- 
plates have a high overlap or fitting factor with the cor- 
rect waveform, which is why the ppE templates have been 
built as generalizations of GR approximations with that 
particular property [60j . The modeling error in the ppE- 
GR templates could be determined by computing fitting 
factors between them and the corresponding GR tem- 
plates, folding in the detector noise curve, or following 
the formalism laid out in This would then provide a 
measure of the efficacy of the ppE templates in detecting 
deviations from GR, or constraining alternative models, 
with the given event. 

Although the practical implementation of an efficient 
post-detection analysis pipeline will be relegated to fu- 
ture work, one should not necessarily restrict such anal- 
ysis to a frequentist approach. An alternative would be 
to perform a Bayesian model selection study |£l|, [62j to 
determine whether the evidence points toward GR or a 
GR deviation (through the non- vanishing of ppE param- 
eters). The idea is that more parameters in a template 
bank will generically improve the fit to the data, even if 
there is no signal present. Thus, any search should be 
guided by careful consideration of priors (such as those 
discussed earlier) or through the computation of evidence 
associated with a given model, for example as in [63| . 
Such considerations are particularly important for low 
SNR events, such as those expected for ground-based de- 
tector sources. 

A clear path then presents itself for possible strategies 
that one might pursue in the study of fundamental bias: 

(i) Given a GR signal and a ppE template, how well can 
the latter extract the former? how much modeling 
error is intrinsic in the ppE template family? 

(ii) Given a non-GR signal and a GR template, how 
much fundamental bias-induced systematic error is 
generated in the estimation of parameters? Can the 
signal even be extracted? 

(iii) Given a non-GR signal and a ppE template, how 
well can the latter extract the former? How well 
can intrinsic and ppE parameters be estimated? 

These questions constitute a starting point for future 
data analysis studies to investigate the issue of funda- 
mental bias in GW astronomy and refinements of the 
ppE framework. 



6 



G. Organization of this paper 

An outline of the remainder of this paper is as foUows: 
Sec. [TT] discusses the anatomy of binary BH coalescence 
in GR; Sec. IIIII derives some model-independent modifi- 
cations to the inspiral GW amplitude and phase, while 
Sec. lIVI does the same for the merger and ringdown phase; 
Sec.|V]uses the results from the preceding sections to con- 
struct an particularly simple example of a ppE template 
bank; Sec. IVII concludes and points to future research. 
Appendix |^ presents a brief summary of well-known al- 
ternative theories and their effect on the GW observable, 
and Appendix |B] describes suggested "exotic" alterna- 
tives to BHs. 

We follow here the conventions and notation of Misner, 
Thorne and Wheeler [4§|. In particular, Greek letters 
stand for spacetime time indices, while Latin letters stand 
for spatial indices only. Commas in index lists stand for 
partial differentiation, while semi-colons stand for covari- 
ant differentiation. The Einstein summation convention 
is assumed unless otherwise specified. The metric sig- 
nature is (—,-(-, -|-, -|-) and we use geometric units where 
G = c = 1. In the Appendices, we reinstate the pow- 
ers of G and c, as one must distinguish the GR coupling 
constants from those in alternative theories of gravity. 



that this classification is somewhat fuzzy, as studies of 
numerical simulation results have shown that a s harp 
transition does not exist between them in GR [s^lssllSq. 
In the remainder of this section, we concentrate on GR 
coalescences and leave discussions about possible modifi- 
cations to the next section. 

In describing the waveform, we shall concentrate on the 
Fourier transform, of the GW response function, 

h{t), as the former is more directly applicable to GW data 
analysis (see eg. [67|, for a data analysis review) . The 
response function of an interferometer, h{t), is a time- 
series defined via the projection h{t) = F+h+ + Fxhx, 
where -f+,x are beam-pattern functions that describe the 
intrinsic detector response to a GW, while /i+,x are the 
projections of the GW metric perturbation hij onto a 
plus/cross polarization basis. We will split the Fourier 
transform of the response function h{f) into three com- 
ponents, ^i(/),/im(/) and hjio{f), corresponding to the 
inspiral, merger and ringdown phases respectively. We 
will further decompose each of these complex functions 
into their real amplitude and phase components: 

^i,m,rd(/) ^ ^i,m,rd(/) 6 i,m,rd(/)^ 



II. ANATOMY OF A COMPACT OBJECT 
BINARY COALESCENCE IN GR 

Let us now consider the quasi-circular binary coales- 
cence of non-spinning HH-like compact objects — namely, 
events that exhibit inspiral, merger and ringdown phases 
akin to binary BH mergers in GR. Certain classes of "ex- 
otic" horizon-less compact objects, such as boson stars 
(see Appendix |B]), may exhibit similar merger wave- 
forms, and insofar as they are viable candidates to ex- 
plain compact objects in the universe that today are clas- 
sified as BHs, they are also within the scope of the ppE 
templates we will construct here. Neutron star mergers, 
for example, are not to be consider part of this fam- 
ily, as these are a separate class of GW sources that 
would require a different ppE expansion (though osten- 
sibly the inspiral phase of a NS/NS merger ppE tem- 
plate would be the same as the BH/BH one, and if the 
merger regime for the former happens outside the detec- 
tor's sensitivity window, as is expected to be the case 
with LIGO/GEO/Virgo, then the binary BH ppE tem- 
plates would suffice). 

We divide the coalescence into three stages: (i) the 
inspiral; (ii) the plunge and merger; (iii) the ringdown. 
In the first stage, the objects start widely separated and 
slowly spiral in via GW radiation-reaction. In the sec- 
ond stage, the objects rapidly plunge and merge, roughly 
when the object's separation is somewhere around the 
location of the light-ring. In the third stage, after a 
common apparent horizon has formed, the remnant rings 
down and settles to an equilibrium configuration. Note 



all of which will be modified in the ppE framework. 

The inspiral phase can be modeled very well by the 
restricted PN approximation (see eg. [6^). In this ap- 
proximation, the amplitude of the time-series is assumed 
to vary slowly relative to the phase, allowing its Fourier 
transform to be calculated via the stationary-phase ap- 
proximation (gqI . FtoI [tH . Improving on the restricted ap- 
proximation leads to sub-leading amplitude corrections 
that introduce higher harmonics in the Fourier transform, 
though we shall not consider that here [7l|, It^]. However, 
we should note that if considering extreme-mass ratio in- 
spirals (EMRIs) or eccentric orbits, then higher harmon- 
ics will play a significant role and should be accounted 
for (see eg. [69|). 

In GR, within the restricted PN approximation the 
Fourier transform of the response function can be written 
as 

Aif) = aj-y^ (3) 
«-:(/) = 2nfto + cPo + J2^ru^''-'^^\ (4) 

fc=0 

where we recall that u = irAif is the reduced fre- 
quency parameter, A4 = Mrf/^ is the chirp mass, 77 = 
mi 7712 /M^ is the symmetric mass ratio, M = mi +m2 is 
the total mass and / is twice the Keplerian orbital fre- 
quency. The coefficient a/ depends on the chirp mass, 
the luminosity distance and the orientation of the bi- 
nary relative to the detector (parameterized through the 
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inclination l and polarization angle 6) via^ 



(5) 



where Q{i, 0) = F+ cos26'(l + cos^ t)/2 - Fx sin 29 cos t 
is a function of the beam-pattern functions F+,x , the in- 
clination 6 and the polarization angle 0.^ The quantities 
to and 00 are constants, related to the time and phase 
of coalescence, while ^™ are constant PN phase param- 
eters that depend on the mass ratio. To 1.5 PN order, 
these coefficients are given by 
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1] 7] 



-2/5 



-3/5 



(6) 



and expressions to higher order are also known [ts'I . 

The merger phase cannot be modeled analytically by 
any controlled perturbation scheme. We shall thus treat 
this phase as a transition region that interpolates be- 
tween the inspiral and the ringdown. During this phase, 
numerical simulations have shown that the amplitude 
and the phase can be fit by [b^, [IJ 



AM) = a,f-l'^r^'\ 



(7) 

(8) 



where /im is the frequency of transition between inspi- 
ral and merger and where 0o and tg are new constants, 
and can be related to (pQ and to by demanding continuity 
of the waveform. The frequency dependence of the am- 
plitude in the merger phase, /~^/^,is not universal, but 
specific to the test system (equal mass, non-spinning) 
studied hcrc'^. 

The ringdown phase can be modeled very well through 
BH perturbation theory techniques (see eg. [t^] for a 
recent review). In this approximation, the GW re- 
sponse function can be written as a sum of exponen- 
tially damped sinusoids, the Fourier transform of each 
being a Lorentzian function. Furthermore, for equal mass 
mergers, a single (least damped) QNR mode dominates 



When modeling the Fourier transform of the response function 
for LISA, one must multiply aj by a factor of V3/2 to account 
for the geometry of the experiment. See eg. Eq. (4.18) and (5.6) 
of [69ll . Moreover, the beam pattern functions further depend on 
the motion of the detector around the Sun, which can be modeled 
as in eg. [zllll 

Note that the function Q as defined here differs from that defined 
in |69|I by a factor of 4. 

Though note that in [76j it was suggested that the exponent 
might be closer to —5/6. 



the waveform, and so to good approximation we can 
write [601 



^Rd(/) — ^l/lM ^ /a4RD 

*nD(/) = 00 + 2^/to, 



1 + 4: 



I-h 



(10) 



where a is the width of the Lorentzian that essentially 
describes the mean duration of the ringdown phase (the 
damping time). Notice that during the plunge and the 
ringdown the phase remains the same, modulo non-linear 
effects, which conveniently guarantees continuity of the 
phase across this boundary. 

The phenomenological parameters (im, /mrd, f) can 
be fit to numerical simulations to find 601 
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^ SttM' 


0.05477^ 


+ 0.0907; 


^-0.19 


1 




TT M 




47rA/' 


0.51r;2 4 


-0.077774 


0.022 


1 



TT M 



(11) 



where the arrow stands for the approximate value for 
equal mass binaries. The quantities (/im,/mrd) de- 
note the boundaries of the inspiral-merger region and 
merger-ringdown region, at which the piece-wise func- 
tion has been constructed to be continuous. Notice that 
the merger and ringdown phases are both rather short 
compared to the inspiral phase, lasting approximately 
(87rAf)"i Hz and (147rM)"^ Hz respectively. Due to the 
short duration of these phases, they can be modeled via 
the simple linear relation in the phase shown above, and 
in fact, as the merger phase is fairly constant, it could be 
absorbed in the overall amplitude. 

The parameterization presented here is purposely sim- 
ilar to that proposed in In that study, however, 
the phase was modeled in all stages of coalescence by the 
same function, namely Eq. @ , where the parameters ■0™ 
are not fixed to the PN value, but instead fitted to numer- 
ical relativity waveforms (see eg. Table 2 in [1^). If one 
prefers, one could model the GW phase in the merger and 
ringdown with the fitted parameters presented in Table 2 
of [60] • Alternatively, one could also follow the prescrip- 
tion of [soj and construct the Fourier transform of the 
response function for an infinite sum of modes, absorb- 
ing the phase into a complex amplitude and paying more 
careful attention to the response function, but we leave 
such details for future work. 



III. MODEL-INDEPENDENT MODIFICATIONS 
TO THE INSPIRAL PHASE 

We shall now consider leading- order deviations from 
GR in the inspiral phase. In the restricted PN and sta- 
tionary phase approximations, the Fourier transform of 
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the response function can be written as (see eg. Eq. (4.5) 
in [H) 

W) ^ A(/)e^*^(/) = rlM 1 e**(//2), (12) 



where A(to) is the time-domain amplitude evaluated at 
the stationary point, defined via ^{to) ~ vr/, and over- 
head dots stand for partial differentiation with respect 
to time. F is the orbital frequency, and F is its rate of 
change evaluated at the stationary point F{to) — //2; 
recall that / is the GW frequency. The quantity ^'(i^) 
is the GW phase, which is given by (see eg. Eq. (4.8) 
in [6|) 



2n 



//2 / f 

dF' {2- — 
F' 



-{F'), (13) 



where t{F) = F/F. The structure of Eq. (fT^ is generic 
for any Fourier transform with a generalized Fourier in- 
tegral that contains a stationary point [70j . 

The orbital frequency and its rate of change completely 
determine the Fourier transform of the response func- 
tion in the stationary phase approximation. The rate of 
change of the orbital frequency is generically computed 
through the product rule F = E {dE/dF)~^ , where E{F) 
is the orbital binding energy. The quantity E is calcu- 
lated by invoking the so-called balance law. This law 
states that the amount of energy carried away in GWs is 
equal to the amount of orbital binding energy lost by the 
binary system, i.e. E = — £gw, where £gw is the GW 
luminosity. The luminosity can be computed perturba- 
tively in terms of a multipolar decomposition of the effec- 
tive energy-momentum tensor of the system (ts} , which 
to leading order yields E ^ \d^Iij/dt^\'^, where lij is the 
reduced quadrupole moment tensor of the binary. 

We will consider GR modifications of the Fourier trans- 
form of the response function during the inspiral arising 
either from changes to the conserved binary binding en- 
ergy (the Hamiltonian, see eg. (79j for a more precise 
definition), or changes to the energy balance law. The 
former are essentially modifications to the conservative 
dynamics, while the latter modifies the dissipative dy- 
namics. We shall tackle each of these corrections next. 



A. Modifications to the Binary Hamiltonian 

The Hamiltonian for a binary system, or equivalently, 
the binary's center of mass binding energy, can usually 
be expressed as a linear combination of kinetic energy 
T and potential energy V: E = T{p) + V{q,p), where 
q and p stand for generalized coordinates and conjugate 
momenta. For a binary system, one can extract the New- 
tonian contribution and rewrite such a generic Hamilto- 
nian as 



E 



V 

y 



M 



ri2 



(14) 



where w = juj — fjl is the relative velocity, fi = ijAI is the 
reduced mass, ri2 is the orbital separation, and M and 
T] are as before the total mass and symmetric mass ratio. 
The quantity V{v,ri2) is some additional contribution 
to the Hamiltonian that can be due either to relativistic 
effects, or to non-GR contributions, for example from ad- 
ditional scalar or dipolar degrees of freedom. Of course, 
not all functions V are allowed, since binary pulsar ob- 
servations have already verified GR in a quasi- weak field, 
which translates into the requirement V ^ MjrYi ^ 1 
when w <C 1. 

From such a Hamiltonian, one can derive a modified 
version of Kepler's third law (or equivalently Newton's 
second law), namely 



M 



1 



ri2 



^{v,ri2) 



(15) 



where lu is the binary's orbital angular frequency and 
il = dV/dr 12. We see then that modifications to the 
Hamiltonian and to Kepler's law are not independent, as 
expected. 

Kepler's law must be inverted to obtain the separation 
as a function of angular frequency, which can then be in- 
serted into the Hamiltonian to obtain dE/dF. Although 
the unspecified form of f2 does not allow an explicit in- 
version, we can always write the inversion implicitly as 



ri2 



(16) 



where x is some function of lo determined from the inver- 
sion of Eq. P3|). As for {V, fl), the quantity x must also 
satisfy the condition x ^ M^/'^w"^/^ in the weak field. 

One possible parameterization of x that satisfies the 
above condition is 



= [1 + Ai xP] 



(17) 



where x = (Mu))^^^ is a PN parameter, while Ai and p are 
dimensionless, real numbers. The parameters Ai and p 
determine roughly where the non-Newtonian correction 
to Kepler's law begins to become important, and they 
must be real since uj is real. Furthermore, the condition 
that X M^/^cj"^/^ when a; ^ 1 imphes that p > in 
the inspiral phase. The choice of parameterization of x 
in Eq. pT|) is equivalent to choosing a parameterization 
for the potential V in Eq. ; with our choice one finds 
that 



V = - 



3Ai 
1 



(18) 



to first order in x. 

We have here chosen a parameterization for x that 
allows an inversion analytically in terms of known func- 
tions. A more general parameterization would be to add 
N terms with exponents di^N > rfz,w-i > ... > p > 0. 
However, the higher the exponent, the smaller its contri- 
bution, since the inspiral phase must terminate outside 
the light-ring where in GR x < 0.4 (for simplicity, we 
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here work with one exponent only). In fact, in GR one 
finds that perturbatively 

1 



ri2 = Mx~U1 + -x{t]-3) 



24257 



37 2 

2520 ^ ^ 12^ 



2 3 41 2 22 ,/a:o\ 

81^^ +192^" +y^'^n-J 



(19) 



when X 1 from Eq. (190) of where here xq is a 
gauge parameter associated with harmonic coordinates. 
From the GR expression, we also infer that the constant 
Ai could depend on the system parameters, such as the 
mass ratio. Also, Ai must depend on some non-GR cou- 
pling intrinsic to the alternative theory, such that in the 
limit as this parameter vanishes, one recovers the GR 
prediction. 

The binary's center-of-mass energy can now be ex- 
pressed entirely in terms of the orbital frequency. Re- 
quiring the circular orbit condition v = ri2Uj, we find 



Ix [1 + XixP]^ - X [1 + XixP]'^ + V{x). (20) 



In principle, we should not expand in a; <C 1, since we 
are interested in corrections that could arise in the late 
inspiral (though we know that even close to merger and 
in GR X < 0.4). To simplify the following analysis, we 
shall perform this expansion regardless, as it only serves 
to motivate the final ppE parameterizations, where x can 
take large values in the late inspiral. Performing the 
X <ti I expansion, we find 



E _ 1 
/7 ^ 2^ 



1 



6 



l+p 



- 4 XixP 



(21) 



and differentiating with respect to F = lu/ (27r) we obtain 
dE 2tt 



dF 



-i^M^x-^l^ [l + 2Ai(l-2p) xP] 



(22) 



The above considerations provide motivation to ex- 
trapolate a model- independent, non-GR parameteriza- 
tion for dE/dF: 



(23) 



where we recall that u = irMf = 2'n:M.F = x^l'^-rf'l^ for 
circular orbits. We have simplified the above parameter- 
ization by introducing A2 = 2Ai(l — 2p) and q = 2p/3. 
Such an expression resembles the PN expansion of GR, 
which one can derive from Eq. (194) of [79( 1 : 



(FE_ 
dF 



27r 
"3 



6 



+ ^4/3^-4/5 



81 57 1 2 



Similarly, the parameter A2 could depend on 77 but it 
must also be proportional to a positive power of the al- 
ternative theory couplings, such that when these vanish, 
one recovers the PN expansion of the GR expectation. 



B. Modifications to the Energy Balance Law 

The energy flux at infinity is usually calculated in GR 
by assuming that the only source of radiation are GWs, 
so that E = —Ccw, where we recall that Cqv^ is the GW 
luminosity. In alternative theories of gravity, however, 
there could be additional sources of radiation from other 
propagating degrees of freedom. The balance law then 
becomes E = -£gw - -Cothor- 

The calculation of the GW luminosity is generically 
performed in two steps. First, one employs perturba- 
tion theory to determine the effective (Isaacson) GW 
stress-energy tensor [sol Isi! ]. the time-time component 
of which provides an expression for £gw in terms of time 
derivatives of the metric perturbation hij. In GR, this 

Second, one 
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expression is schematically Cc\ 
performs a multipolar decomposition of the metric into 
mass and current multipoles '78l|, which are given by in- 
tegrals over the sources and for point sources are simply 
proportional to their trajectories and velocities. In GR, 
one then schematically finds hij ~ lij / R, where is the 
reduced quadrupole moment and R is the distance to the 
observer. Combining both steps into one single proce- 
dure, one finds that the GW luminosity in GR is propor- 
tional to third time-derivatives of the reduced quadrupole 
moment Low oc I ij I^'' ■ 

In general, one can decompose the GW luminosity 
in a superposition of mass and current multipole mo- 
ments 17811 



Cgw « kmM^ + kdD,D' + Kqiij r 



(25) 



where Ai is the monopole of the system, D is its dipole 
moment and km,d,q are proportionality constant for the 
monopole, dipole and quadrupole contributions respec- 
tively. The quantity km has the effect of re-normalizing 
Newton's gravitational constant G, and thus, can be bun- 
dled together with kq. On the other hand, the quantity 
kd is related to the self-binding energy of the gravitat- 
ing bodies, and it produces an independent contribution 
to the GW luminosity. In GR, these two contributions 
(monopole and dipole) identically vanish due to center- 
of-mass and linear-momentum conservation. Evaluating 
Eq. (j25p with standard definitions of multipole moments 
for a binary system 78], one finds to leading order 



E : 
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15 



(KlV 



-kdG 



other 5 



(26) 

where ^1^2 are the so-called Peters-Mathews parame- 
ters, whose numerical value depends on the alternative 
theory under consideration, while kjj is a dipole self- 
gravitational energy contribution, with Q the difference 
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in self-gravitational binding energy per unit mass (see 
eg. Eq.(10.84) in |/7,]). For example, in GR {ki, K2, kd) — 
(12, 11,0), while in Brans-Dicke theory without a cosmo- 
logical constant K2, hd] = [12 — 5 (2 + ujbd)^^, 11 — 
11.25 (2 + W3d)-\2 (2 + Wbd)-'] 0, where cJbd is the 
BD coupling parameter. 

Using the modified Kepler's law from the previous sub- 
section, we can re-express the power as a function of fre- 
quency only, namely 



5 



96 



Co 



the 



(27) 



where we have set r — Q, corresponding to a circular or- 
bit. The dipolar term dominates over GR's quadrupolar 
emission during the early inspiral, which explains why 
Brans-Dicke theory, whose dominant signature is pre- 
cisely dipolar radiation (see Appendix IA 1[) , has been con- 
strained so well by binary pulsar observations. On the 
other hand, the modification to Kepler's law becomes im- 
portant as the binary approaches the end of the inspiral 
phase, as expected since it is meant to model dynamical 
strong-field corrections to GR. 

The quantity £othcr serves to parameterize both ad- 
ditional contributions from scalar or vectorial degrees of 
freedom, as well as high-order curvature corrections to 
the GW effective stress-energy tensor. In the case of the 
former, one could for example have £othcr — C^p oc 0^, 
where is some scalar field that is activated when the 
binary approaches the merger phase. In fact, such is 
the case in dynamical CS modified gravity (see Ap- 
pendix IXS)) . where the CS field becomes strongly sourced 
during the late inspiral, plunge and merger, as the Rie- 
mann curvature grows. Here, one could for example have 
-Cothnr = Cr*r oc • (Vfe/;„i)^p, an interaction that 



-'Other 

is absent in GR. 



Modified GW Phase 



We can now combine the results from the previous two 
sections to find F: 



F 
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'10/3 



(28) 



where we have rescaled ki = kd = ^i^dG V 
C 



;2„2/5 



other 
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32 



'Cothcr, and where the new constant parameter 
A3 can be related to (ki , Ai , p, 77) , though the expression is 
not particularly illuminating. We have also here dropped 
second-order terms, such as the influence of the modified 
Keplerian relation on dipolar emission. The parameter Ri 
has the effect of re-normalizing the value of G or c, which 
would not be directly observable in a GW detection, as it 
would be directly reabsorbed by physical quantities, like 
the chirp mass (unless one possesses an electromagnetic 
counterpart). The parameter Rd parameterizes dipolar 



deviations and depends on the sensitivities of the bodies. 
The parameters ki and q parameterize deviations from 
Kepler's law, while £othcr describes possible additional 
sinks of energy. 

We can now integrate Eq. (fT3)) to obtain the modified 
GW phase: 



-- ~(/)c + 27r/ie 
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7) {u')-^^'^Coth..dF', 



where A4 and Rd are new parameters that again can be 
related to the old ones, though these relations are unil- 
luminating. As is customary, we have here linearized in 
the GR corrections to be able to perform the integral 
analytically, an approximation that is valid in the inspi- 
ral, although not necessarily so in the plunge and merger. 
Note that the term proportional to Rd mimics the Brans- 
Dicke correction to the GW phase, which is known to 
scale as m^^/"^ (see discussion around Eq. (IA6|) ). while 
q = —1 mimics the effects of a massive graviton on the 
propagation of GWs, which is known to scale as (see 
discussion around Eq. (jAlOp V 



D. Modified GW Amplitude 

The GW amplitude depends both on the amplitude 
of the time-domain waveform evaluated at the station- 
ary point and the rate of change of the frequency F . 
The latter has already been parameterized in Eq. (|28p . 
The former, however, requires further study. In GR, the 
quadrupole formula requires that 



(30) 



which leads to the following waveform amplitude to lead- 
ing order for a binary system 



(31) 



where Q(l, /3) is as defined after Eq. ([5]). Evaluating this 
amplitude at the stationary point and using the modified 
Keplerian relations of Eq. (|17p . we find 



AM 
Dl 



(32) 



Combining this expression with Eq. (|28p . we find that 
= is equal to 



Qii,P)f-'/' (Ri 
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where ki is another constant and we now kept only the 
leading-order corrections. Such an expression agrees 
identically with the GR expectation when Ri — I and 
all other modified parameters vanish. 

The above results hinge on the assumption that the 
quadrupole formula in Eq. ((30)) holds to leading order 
in the GR modification. In well-known alternative the- 
ories of gravity, including scalar-tensor ones, Eq. ([50]) 
is indeed valid to leading order, acquiring subdominant 
corrections. These corrections can be modeled by the 
undetermined function Mother, in the same fashion as in 
Sec. IIIIBI Thus, the generalization of Eq. ([50]) leads to 
redundancies in the parameterization that we are after. 
In Sec.|V]we shall use the insight gained from the analysis 
of this section to infer a generalized parameterization for 
a model-independent modification of the inspiral gravi- 
tational waveform. 



IV. MODEL-INDEPENDENT MODIFICATIONS 
TO THE MERGER AND RINGDOWN PHASE 

A. Merger Phase 

Let us now consider leading- order deviations from GR 
first during the merger and then the ringdown phase. The 
merger phase is the least understood regime of a binary 
coalescence, because no controlled approximation scheme 
can be employed to solve the field equations. In GR, 
numerical simulations of the merger of two, non-spinning, 
equal-mass BHs has produced the phenomenological fit 

hUf)ocu-^^'e'^-^f\ (34) 

where the phase is a linear function '^A-iif) = 6 + eu. 
Comparing to the results of Sec. [ll] we find that 5 = (j)o 
and e = 2tQA4^^. Linearity in the phase can be assumed 
due to the short duration of the merger phase (usually 
less than 2 GW cycles for BHs). This functional form is 
an approximate fit and found only a posteriori, i.e. after 
a full numerical solution has been obtained. 

Lack of formally derivable analytic approximations of 
the merger waveform is not exclusive to GR. Further- 
more, due to the dearth of numerical studies of this 
regime in alternative theories (but see [13] for promising 
plans to explore mergers in scalar- tensor theories), little 
guidance can be found on how the merger regime may 
differ from GR. How do we then parameterize model- 
independent corrections of the merger waveform when 
we possess no analytic or numerical guidance? For lack 
of a better prescription, we will here assume that the 
waveform produced during the merger can adequately be 
described as an interpolation between the inspiral and 
the ringdown waveforms. This is a good approximation 
to GR waveforms, and hence is a rather conservative ap- 
proach to a ppE extrapolation. Specifically then, we pro- 



pose 

oc 

/iM(/)-^7M,„u^"-e^^*+^"\ (35) 

where 7m, « and CM,n are an infinite set of undetermined 
parameters. We shall here keep only one term in the poly- 
nomial, setting 7m, = 1, Cm,o = c and all other amplitude 
parameters to zero, while we shall leave e undetermined 
in the phase. The parameters (7, S) are set be requir- 
ing continuity between the inspiral-merger interface at 
u = UjM- 7 = A(/im)w™ and 5 = *i(/im) - cUim. The 
choice of a single term in the power-series representation 
of the Fourier transform of the merger response function 
forces us to consider coalescences with prompt mergers. 
Delayed mergers can occur, for example, with exotic com- 
pact objects such as boson stars [83|, in a process akin 
to delayed-collapse scenarios in NS mergers. Of course, 
additional terms could be added to the series of Eq. ([55)1 
to account for such scenarios, however for simplicity here 
we will restrict attention to prompt mergers. 

B. Ringdown Phase 

Let us now consider the leading- order modification to 
GWs that could be produced during the post-merger 
phase due to alternative theory corrections. GWs emit- 
ted in this stage are produced by oscillations in the rem- 
nant spacetime as it settles down to its final equilibrium 
configuration. Simulations show that in GR the ringdown 
very quickly becomes linear, and we will assume the same 
here for the ppE expansion. In GR, the metric perturba- 
tion in general is a sum of quasi-normal ringdown (QNR) 
modes, which are exponentially damped sinusoids, and 
so-called tail terms, which are constant phase solutions 
that decay via a temporal power law. The QNR modes 
dominate the early-time post-merger waveform, which is 
why we label this phase the ringdown phase. 

We can write the ringdown time-domain response func- 
tion as 

K„{t) = AaD(t)cos[^'Ho(i)]: (36) 

where the phase is purely real and the amplitude is a 
decaying function of time. Due to the assumed linearity, 
we can then decompose the amplitude and phase into 
a sum of harmonics h^^ = ^rd^IO cos[*S^™(0], 

where (£, m) labels a spheroidal harmonic mode, and n 
is an integer describing the overtone of the given mode. 
For the QNR component of each mode, the amplitude is 
a decaying exponential 

QNR^n^m^^^ = 6'(t) A^^^e-*/^"^" , (37) 

where Tnim is the decay constant of the mode, 6'(-) is 
the Heaviside function that guarantees the waveform is 
well defined at negative infinite time, and AJf^™ is some 
constant amplitude determined by continuity with the 
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merger phase, which as such must depend on the beam- 
pattern functions. The phase of each QNR mode is 

= 27rffrnnt, (38) 

where fimn is the real quasinormal frequency of the 
(n, i, m) mode. 

In GR, the Fourier transform of the QNR component 
of the ringdown waveform is thus a sum of Lorentzian-hke 
functions: 

nl-m 



This result can be obtained by rewriting the cosine in 
Eq. (j36|) as a sum of complex exponentials using Euler's 
formula, and then following the doubling prescription 
of 84], where one essentially replaces t/Timn \t\/Temn 
in the Fourier integral and then multiplies the result by 
a correction factor that accounts for the doubling [30| . 

Tail terms in the post-merger waveform are induced 
either by linear interactions between the GWs and the 
background, or non-linear behavior right after merger. 
Linear tails have a long history in GR [s^ [s^, [s^l and 
they are controlled by the asymptotic form of an effec- 
tive potential far from the source. Such tails depend 
on the harmonic one is looking at but generically for a 
Schwarzschild background they can be parameterized by 

Recently in (s^l it was suggested that "non-linear tails" 
could be produced during merger. Such tails are struc- 
turally similar to the linear ones, except that one replaces 
the power 2£ + 3 b, where b is some positive integer. 
In this case, the Fourier transform is given by 

-"/...(/) = y: — (41) 

for positive frequencies (the Fourier transform for a 
Schwarzschild linear tail can be easily inferred from the 
above equation). In GR then, one expects the complete 
ringdown GW to be a superposition of QNR terms and 
tail terms, leading to a Fourier transform that is the lin- 
ear combination of Eq. p9|) and (|4T|) . 

What generic type of modifications could we expect 
from alternative theories of gravity? One possibility 
would be to modify the power in the exponential damp- 
ing part of the ringdown, i.e. in Eq. (|39p exp(— t/T„fm) 
exp(— t'^/r^^^), with a positive real number. However, 
a ^ 1 does not seem to be very "natural" , and would 
significantly complicate the analysis. For this study we 
simply keep a = 1. Another possible modification is to 
correct the tail behavior of the ringdown waves. Linear 
tails, however, are predominantly controlled by the dom- 
inant 1 /r behavior of an effective potential in the far-field 



field, which in turn depends on the 1/r behavior of the 
background metric tensor. If we require that all stable 
post-merger metric tensors approach the Schwarzschild 
solution at late times in the far-field (by Solar System 
tests), it is unclear how much linear tails would be mod- 
ified by non-linear terms in the gravitational field. On 
the other hand, non-linear tails could be modified in 
alternative theories of gravity, as suggested in ^8Sl, 93] ■ 
Nonetheless, alternative theory corrections to either lin- 
ear or non-linear tail terms can both be parameterized 
by some undetermined power 5, leading to the Fourier 
transform shown in Eq. (PT|) . 

For simplicity then, we will propose a modified ring- 
down waveform with a similar functional form as the GR 
result given in Eq. ((39l) . though allow the decay con- 
stants, frequencies and powers to take on non-GR values. 



V. PARAMETERIZED POST-EINSTEINIAN 
FRAMEWORK 

In this section we collate the results from the previ- 
ous sections to construct ppE extensions of the inspi- 
ral, merger and ringdown waveform segments of a binary 
merger. In each of these phases, we shall divide this task 
into the construction of a ppE phase and a ppE ampli- 
tude, such that the waveform is given by: 

^^r(/)exp[zfr(/)], (42) 

where A = (I, M, RD) for the inspiral, merger and ring- 
down phases respectively. Then, we shall put these wave- 
forms together to construct a ppE waveform for the entire 
coalescence. We conclude this section with a discussion 
of the data analysis cost incurred due to the addition of 
ppE parameters to the waveform. 



A. Inspiral 

Let us employ the insight gained by the analysis of 
Sec. mil to propose inspiral ppE waveforms. Consider 
then the ppE GW phase, a generic form of which can be 
extrapolated from Eq. (|29[) . One simple parameterization 
for the test system is the following: 

7 

*r (/) ^-\-c^c+ 2nfto + pu" + u-'/^ Yl ^r^'/', 

(43) 

where ip™ are PN phase coefficients, some of which are 
given in Eq. ([S]), while (3 and b are ppE parameters. Such 
a parameterization allows one to easily map between BD 
GW modifications and massive graviton ones, through 
the parameter choices (/9, 6) — (/3bd,— 7/3) and (/3, 6) — 
(/3mg,-1) respectively. 

One generalization would be to allow more than one 
alternative theory to be modeled simultaneously by the 
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ppE waveform. With the test system in mind, such an 
idea suggests 



(44) 



The benefits of this parameterization is that it exphc- 
itly incorporates the GW phases predicted by alternative 
theories, as well as the PN prediction, while allowing, 
through the %^ term, for new alternative theory modi- 
fications that have not yet been developed. Note how- 
ever that now we have introduced a degeneracy into the 
waveforms; i.e. when h = —T/i, waveforms along the line 
7 -I- /3 = const, are identical. In a sense then this is a 
non-optimal expansion. 

A better suggestion might then be to postulate 



ppE,3 



(/) 



TT 



2^ ft. 



-7/3 



Y^ct^ku''"'. (45) 



fe=0 



The controlling factor in the phase expansion has been 
modified to —7/3 instead of —5/3 to allow for dipole radi- 
ation. The phase parameters (pk can take on any values; 
the ppE-GR templates will be those where (t>k = V'™- 
Continuing in this vain then, the most general expansion 
we propose here for the phase is the following: 



4- 



77 
"I 



N-1 



2t: ft, + J2 



(46) 



fe=o 



PN 

k ' 



where now i^^ and bk are each a set of N ppE parame 
ters. The ppE-GR templates are those with (j)k = '0, 
and bk = {k — 5)/3. Again, this family can also capture 
the inspiral phase of several well known alternative the- 
ories (see Appendix A) with appropriate choices of the 
parameters. The parameterizations in xj/^p'^''^'^ are a gen- 



eralization of the studies in [25|, l56|, l57[ , except that we 
here allow for non-GR exponents in the u-dependence of 
the phase that could model dynamical strong-field cor- 
rections to GR. 

Now let us concentrate on amplitude parameteriza- 
tions. A generic form for the amplitude that can be 
extrapolated from Eq. ((33)) is the following 

Arif) = (I) TT-y^^Qf-'^' [1 + aix^onK] 

(47) 

The quantities a(AGR) and a are a ppE function and pa- 
rameter respectively, where the former could depend on 
the GR parameter vector Agr = [M, l, F+^x)]- Such a 
non-GR deformation of the GW amplitude is not de- 
generate with high-order PN corrections, since these 
scale with higher-harmonics of the GW signal, instead of 
the fundamental mode. Moreover, such a parameteriza- 
tion reduces to well-known alternative theory predictions, 



such as CS modified gravity when a = 1 and a reduces 
to the second term in Eq. (|A18p . Of course, one could in 
principle introduce N terms to the GW amplitude with 
different frequency exponents, though for simplicity we 
will not do so here; furthermore, GW detectors are most 
sensitive to phase, so there might be less to gain by in- 
troducing expansive amplitude modifications. 



B. Merger and Ringdown 

Here we employ the considerations discussed in Sec. IIVI 
to construct a ppE waveform that models the merger and 
ringdown phase in a non-GR, model-independent fash- 
ion. For the merger phase, we have (see Eq. (j35p with 
Eq. gSD) 



(48) 



The ppE parameters during the merger are then the real 
numbers (c, e), while the parameters (7, S) are set by re- 
quiring continuity with the inspiral phase. When the ppE 
parameters are (c, e) = (—2/3, 2ioAl~^) one recovers the 
GR prediction for the test system. 

A simple proposal for a ppE parameterization of the 
ringdown waveform is 



l+4TT^T^K{f- fno) 



(49) 



with ^'gjf = 0, and where /hd and r are the dominant 
quasi-normal (QN) frequency and decay time, k is a real 
ppE function of the final mass of the remnant, d is a real 
ppE parameter, the quantity C is determined by continu- 
ity with the merger phase, and where we have chosen the 
transition to occur at the fundamental ringdown mode 
/mrd = /rd- This models the ringdown with a single 
decaying mode, although even in GR there are an infi- 
nite number of such modes excited during ringdown (a 
possibility we analyze below). 

A more general ringdown waveform is as follows (oij 



■if) 



c- 



Tlr, 



1 



(50) 



where femn and Timn are the frequency and damping time 
of the (n^TO)-th mode, while the ppE function k and pa- 
rameter d are assumed mode-independent. One could in 
principle postulate different ppE parameters for different 
modes, but this would introduce an immense number of 
such parameters, not all of which would necessarily be 
independent. 

Traditionally, suggested tests of alternative theories of 
gravity using ringdown GWs are based on the "no-hair" 
theorem of GR. This theorem states that given the mass 
and angular momentum of a perturbed BH, all frequen- 
cies and damping times of the waveform are determined. 
For example, in GR the (220) mode is given by [t^, Iq^ 



MlU2 

Q 



220 



1.5251 - 1.1568(1 -a) 
0.700+ 1.4187(1 - a) 



0.1292 



0.4990 



(51) 
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where the total mass of the remnant is M, the spin Kerr 
parameter of the remnant is a = J/M^, and we have 
defined the quaUty factor Qimn = ^imn/ {'^Timn)- The 
standard no-hair test assumes that two modes have been 
detected, the first of which can be used to determine a 
and M, while the second one can be used to check for 
consistency through equations similar to those above for 
additional modes. The violation of the two-hair theorem 
is common in alternative models for compact objects (3l| . 
With the parameterization of Eq. (jSOp . no assumption 
has been made about the relation between different tmn 
modes, and thus, such a test is still possible. 

One fact that should be kept in mind is that the 
strength of decay of the QNR modes (i.e. the scaling of 
Timn with harmonic number) depends strongly on how 
"exotic" the compact object is. In GR, the first quasi- 
normal mode is dominant because r^mn increases rapidly 
with and m. This is not the case, for example, with 
boson stars [3l| ; in such objects there are classes of modes 
with comparable damping times, yielding a signal that 
cannot be extracted with one (or a small set) of modes. 
We shall not concern ourselves with this issue here, be- 
yond pointing out that some strategies have been devel- 
oped to deal with it [soj . 



C. The Complete ppE Waveform 

In summary then, the complete ppE waveform we have 
proposed takes on the following piece-wise continuous 
form 

r/ir'^(/) for/</™, 

Hf) = { krif) for < / < /m,o, (52) 

where A is any of A = [1, . . . , 4] and B is any of i? = [1, 2]. 
All pieces of the waveform are decomposed into an ampli- 
tude and a phase with the structure of Eq. (|42| . During 
the inspiral, the amplitude is given by Eq. (|47p . while the 
phase can be any of Eq. (gSl), (gll), (jlSl), (gSl), depending 
on the level of bias one wishes to allow for. During the 
merger, the amplitude and phase are given by Eq. (j48p . 
while in the ringdown the phase is zero and the amplitude 
is given by either Eq. (|^ or ([50)1 depending on the level 
of generality desired. The quantities /im and /mrd are 
the inspiral-merger and the merger-ringdown frequency 
transitions respectively. These quantities can also be con- 
sidered ppE parameters, or alternatively /im and /mrd 
can be set to the light-ring frequency and dominant ring- 
down mode of the corresponding remnant ppE-GR BH, 
for example. 

Note that these waveforms are not analytic; they are 
continuous but non-differentiable at the patching fre- 
quencies between the phases. Of course, if desired, ppE 
families with higher levels of differentiability could be 
constructed, though our main goal here was to introduce 
the ppE framework with a simple (yet potentially useful 
and relevant) model. 



Furthermore, one could easily extend these waveforms 
with additional ppE parameters to try to capture larger 
classes of deviations from GR. However, there is an in- 
creasing price one must pay the more parameters are 
added to a fitting routine: an increase in the number 
of filters needed and a deterioration in the accuracy of 
parameter estimation. The extra computational cost 
needed to deal with additional parameters could be pro- 
hibitively high. For example, one can show that the 
introduction of one single extra parameter to ringdown 
waveforms increases the number of required templates to 
cover the parameter space from 500 to 10^ in a frequen- 
tist, matched-filtering scheme [93i]. In view of this, one 
must ensure that additional parameter included are truly 
well-motivated, as it is explained in the Appendices. 

At first glance, another potential problem with these 
ppE templates, were they to be used directly for GW 
searches of equal mass, binary BH inspirals, is they do 
contain some modeling error, i.e. with all ppE parameters 
set to their GR values, the waveforms will not exactly de- 
scribe the merger event in GR. In principle one can con- 
struct ppE waveforms with as small a modeling bias as 
desired, though this may not be practical, in particular if 
hybrid waveforms using numerical solutions for the late 
inspiral and merger phase of the waveforms are needed. 
The numerical waveform segments will only be known as 
sets of numbers — either direct samplings of the simula- 
tion waveforms, or coefficients of some truncated series 
expansion fit to the latter — and it is unclear how such 
waveforms could then be extended in any well-motivated, 
model independent manner as done here. 

However, as discussed in the Introduction, we do not 
propose to use the ppE-templates for direct detection, 
but rather in post-analysis after a detection has been 
made with GR templates. Some modeling error is then 
not much of an issue, other than possibly limiting the ac- 
curacy with which the given ppE family can measure de- 
viations from GR. Also, there are alternative techniques 
to "brute-force" matched filtering that would alleviate 
some of the problems with increased numbers of parame- 
ters, though we leave further study of this in the context 
of the ppE framework to future investigations 9j| . 



VI. CONCLUSIONS 

In this paper we have considered the notion of funda- 
mental bias in GW astrophysics. The main source of this 
bias is the assumption that GR correctly describes GW 
physics throughout the generation of the waves, and their 
subsequent propagation. If GR is wrong in the dynami- 
cal strong-field regime, a region yet to be constrained by 
observational tests, and templates based solely on GR 
are used, such a bias could cause wrong inferences to be 
made about the population and nature of corresponding 
sources in the universe. 

We have introduced the parameterized post- 
Einsteinian, or ppE framework to help remedy this 
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problem of fundamental bias for classes of sources that 
do not deviate too radically from GR to prevent detec- 
tion with GR templates. The ppE framework introduces 
non-GR parameters in the waveform to model non-GR 
deviations, in analogy to the ppN framework developed 
for Solar System tests. As an example, we focused 
on one source for which there is good observational 
evidence that GR correctly describes (most of) the 
event: the quasi-circular inspiral, merger and ringdown 
of equal-mass, non-spinning BH-like objects. Each 
ppE parameter is ostensibly selected to characterize a 
violation of some fundamental GR property or principle, 
for example emission of dipole radiation; we were 
furthermore careful to choose parameters so that the 
ppE template family would be sensitive to known effects 
of certain theories of gravity, as reviewed in Appendix 
A. 

The acknowledgement of bias and the formal introduc- 
tion of the ppE framework opens up an entire new set of 
possible studies, as discussed at the end of Sec. II Fl De- 
pending on whether the signal is a GR or non-GR one 
and whether the template is a GR or ppE one, one can 
attempt to answer several important questions that are 
critical to data analysis. For example, with a non-GR 
signal and a GR template, one can determine how much 
systematic error induced by fundamental bias might con- 
taminate parameter estimation and perhaps even signal 
extraction (depending on SNR). Alternatively, given a 
non-GR signal and a ppE template one can attempt 
to determine how well GW observations can truly con- 
strain generic GR deviations. Through a systematic re- 
finement of the ppE framework and the development of 
different ppE template families appropriate for different 
sources (spinning BH binary coalescences, extreme-mass 
and intermediate-mass ratio inspirals, etc.), we shall thus 
be able to determine the true potential of GW observa- 
tions to constrain our standard model of the gravitational 
interaction. 

There are numerous additional avenues to be explored 
here. Several include studying the systematic errors 
fundamental bias could introduce on inferences about a 
population of detections made with pure GR templates 
{i.e. question (1) in Sec. II C[) . designing ppE families for 
additional events (larger classes of binary BH mergers, 
binary NS or NS/BH mergers, NS oscillations, etc.), re- 
fining the proposed ppE detection pipeline discussed in 
the Introduction,, and exploring how strongly particular 
alternative theories may be constrained using ppE tem- 
plates. 

In addition, it would be very useful to understand how 
the merger phase in compact object coalescence differs in 
alternative theories. Such an understanding can only be 
gained via numerical simulations of such merger events in 
alternative theories that possess well-posed initial bound- 
ary value problems. This would give some indication of 
whether a simple interpolation between inspiral and ring- 
down would suffice, or if not, give guidance on how to 
model this phase in a ppE fashion. 



The proposed ppE framework could also become more 
powerful if coupled to coincident electromagnetic detec- 
tions, for those events that also emit photons. In such a 
case, some of the system parameters, for example the 
mass and the luminosity distance, could be indepen- 
dently (electromagnetically) determined, thus reducing 
the full (system plus ppE) parameter space and allow- 
ing for a more accurate measurement of ppE parameters. 
Electromagnetic observations, however, might also be af- 
fected by some of the GR modifications discussed here. 
A fruitful direction for future research would then be to 
map the ppE parameters introduced here to modifica- 
tions of electromagnetic observations. 
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APPENDIX A: A BESTIARY OF ALTERNATIVE 
THEORIES OF GRAVITY 

In this appendix we shall provide a brief review of some 
well-known alternative theories of gravity. We refer the 
reader to the extensive references presented in this Ap- 
pendix for further reading. 

1. Scalar- Tensor Theories and 
Jordan-Brans-Dicke-Fierz Theory 

Scalar-tensor theories modify the action by the inclu- 
sion of a scalar field that couples to the Einstein-Hilbert 
term with a field-dependent coupling function A{lp) and 
a potential V{ip), namely (see eg. [a0| for details): 

+ S^U^^^.,A^^)g^,], (Al) 

where commas stand for partial differentiation, g is the 
determinant of the Einstein metric g^^, S'mat is the ac- 
tion for ■(/'„„t7 and we have neglected the cosmological 
constant term. The scalar-tensor theory action as writ- 
ten in Eq. (lAip is sometimes said to be in the Einstein 
frame [95|, i96|] . 

The action as written in Eq. (|Aip might seem to define 
a non-metric theory, because the external matter degrees 
of freedom ■f/'mat couple both to the metric tensor and 
to the scalar field <p [g]. As we explain below, however, 
scalar-tensor theories are indeed metric as the scalar field 
acts only as a mediator that influences the metric, but it 
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is the latter that determines the motion of external mat- 
ter fields. By performing the conformal transformation 
A'^{(p)gfj_i, we can map the above action into 



1 



IGttG 



4>R 



^(0) 



'■V 
(A2) 



which is clearly a metric theory. In Eq. (|A2[) . the new 
field (/) is defined via (j) = A^'^ , the coupling field 
is defined as lu = (a^^ — 3)/2, where a = A^^/A. 
Equation (IA2|) is sometimes said to be in the Jordan 
frame |95|,196|]. 

Variation of the action with respect to the metric and 
the scalar field leads to the modified field equations [7, 
19711 ■ which in the Jordan frame are 



□0 



1 



3 + 2^(0) 
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la 



(A3) 



where □ = g^'^d^n is the D'Alembertian operator and 
T^"^"' is the matter stress-energy tensor. Because of the 
form of the modified field equations, scalar-tensor the- 
ories are sometimes thought of as modifying Newton's 
gravitational constant via G ^ 1/0. When the coupling 
uj{(l)) — cjbd is constant, then the scalar-tensor theory 
reduces to BD theory [98| . 

Scalar tensor theories satisfy all of the necessary cri- 
teria we discussed in the Introduction. First, scalar- 
tensor theories can be written as purely metric theories, 
as shown above. Second, these theories reduce exactly 
to GR in the limit cJbd — > oo, and thus, it passes all 
precision tests for a sufficiently large coupling. As for 
the dynamical strong-field behavior, one might expect 
the gradients of the scalar field to become large close to 
merger, although scalar-tensor theories are not built to 
introduce high-order curvature corrections. Early studies 
of scalar-tensor theories in NSs suggest that the dynam- 
ical strong-field behavior can be greatly modified, a, pro- 
cess sometimes called spontaneous scalarization [22, l99| . 

The additional theoretical criteria discussed in the In- 
troduction are also satisfied by scalar-tensor theories. 
Such theories have be en shown to be well-posed as an 
initial value problem [lOOl Il0l| |. Moreover, BD the- 
ory has been shown to possess stable Schwarzschild- 
like BH solutions that pass all precision tests [l^l for 
sufficiently large coupling. Finally, scalar-tensor the- 
ories are well-motivated from the low-energy effective 
theory limit of string theory. More precisely, when 
one integrates out all the string quantum fluctuations, 
one finds that the higher-dimensional string theoreti- 
cal action reduces t o a local field theory similar to a 
scalar-tensor theory jl02l . Il03j |. The mapping between 
scalar-tensor theories and string theory is (/> = e"^"^, 



where ip is the dilaton field that couples to matter via 
m(V') = m exp(/3'0), where m is a constant and (3 me- 
diates the dilatonic coupling. By requiring the Strong 
Equivalence Principle to hold, one can show that the ac- 
tion takes the form of E g. (|A2I) in the Jordan frame with 
UJ = (I//32 _ i2)/8 [l02[ll03j . 

Scalar-tensor theories can also be mapped to the re- 
cently popular f{R) class of theories. In the latter, 
one modifies the Einstein-Hilbert action by replacing the 
Ricci scalar by some functional of this quantity. Of 
course, such an action defines an infinite class of theories 
as f{R) is arbitrary. Such arbitrariness removes a certain 
amount of predictability from f{R) theories. Moreover, 
one can show that f{R) theories are equivalent to BD 
theory with = 0, via the m apping (j) = df{R)/dR 
and V{(j)) = Rdf{R)/dR - f{R) ^M, M^- Due to its 
reduced predictability and equivalence with BD theory, 
we shall not consider f{R) theories further in this paper. 

The emission and propagation of GWs in scalar-tensor 
theories is different from that in GR. The main effect 
of BD theory is the emission of monopolar and dipolar 
energy, which modifies the GW flux formula via E = 
i?BD + -Egr: where the d omin ant BD correction is dipolar 



(A4) 



and of the form HJ, \M, M 
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where we recall that M = mi + m2 is the total mass, 
ri2 is the binary orbital separation and rj — mim2 / M'^ is 
the symmetric mass ratio. The quantity 5" is a function 
of the sensitivity of the compact objects, related to the 
self-gravitational binding energy per unit mass, which is 
determined by the equation of state. Such a modiflcation 
arises because BD theory also corrects the compact ob- 
ject's effective mass, causing the location of the center of 
gravitational binding energy to disagree with the center 
of inertial mass. The effective mass now depends on the 
internal structure of the bodies, thus violating the Strong 
Equivalence Principle. Although for pure NS binaries or 
pure BH binaries this effect is either small or identically 
cancels, it does not vanish for mixed (BH-NS) binaries. 

The Fourier transform of the GW metric perturbation 
is greatly affected by the BD modification to the GW flux 
formula, since this propagates into modifications to the 
evolution of the orbital radius and GW frequency. One 
can show that the transform of the response function 
(defined in Sec. |TT1) takes the form [TtL [23L [Toej 

h^n = h^ne"^^°, (A5) 

where /igr is the Fourier transform of the response func- 
tion in GR. The Brans-Dicke phase correction, 0bd = 
—Pbo u~'^^^, where we recall that u = irM-f is the di- 
mensionless frequency parameter, the chirp mass M = 
rf/^M and the Brans-Dicke phase parameter 
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neglecting terms of 0{uj~^). Note that m 
Wbd oo one recovers GR. In the Solar System, Doppler 
tracking of the Cassini spacecraft has led to the {2a) con- 
straint cjbd > 4 X 10^ [58]. Studies have suggested that 
a GW detection could place constraints on Wqd that are 
one or two orders of magnitude larger than the Cassini 
one [13, HE HI Hi. 

GWs in scalar-tensor or BD theory are not restricted 
to having only two polarizations, namely, these theories 
also allow for a breathing mode. GW detectors, however, 
are not directly sensitive to breathing modes due to the 
geometry of the interferometer. In spite of this, GW 
detectors will be indirectly sensitive to these modes, since 
they will carry energy away from the binary and thus 
naturally lead to a modification of the energy balance 
law. The correction to the GW amplitude, and thus the 
GW response function, however, is subdominant si nce i t 
leads to ©(I/cjbd) corrections to the GR prediction flOT^. 

Apart from the inspiral phase, the merger and ring- 
down phases are also modified in scalar-ten sor t heories. 
The ringdown phase has been analyzed in jl08l |. where 
it was found that the QNR frequencies in NS mergers 
are shifted in such theories. The plunge and merger 
regime have not yet been studied, although one expects 
the largest GW modifications to arise in this region. In 
fact, [99] has shown that dynamical strong- field gravity 
effects in scalar-tensor theories can greatly affect NS pa- 
rameters and their oscillations. 



2. Massive Graviton Theories 

Massive graviton theories is the name given to mod- 
els in which the gravitational force is propagated by a 
massive gauge boson, i.e. a graviton with mass rUg ^ 
or Gompton wavelength Xg = h/{mgc) < co. In GR, 
GWs can be thought of as massless gravitons propagat- 
ing at the speed of light. If the gravitons possess a finite 
Gompton wavelength however, GWs will travel at veloci- 
ties less than the speed of light. When we refer to massive 
graviton theories in the main text of this paper, we really 
mean corrections to the propagation of GWs as induced 
by an effective graviton mass. 

On a classical effective level, several theoretical frame- 
works exist wher e G Ws p ropagate at speeds less than 
that of light [109, 113, [lllj. One example is Rosen's bi- 
metric theory ,109 | , where gravitons follow null geodesies 
of some fiducial metric 77^^, while light follows null 
geodesies of some other metric 5^1/ @, [3|- A- similar, more 
recent example is Visser's massive graviton theory [lioj. 
where the graviton is given a mass at the cost of intro- 
ducing a non-dynamical background metric, i.e. a prio r 
geometry. Another recent example is TeVeS theory |lll| , 
where the introduction of a scalar and a vector field lead 
to subluminal GW propagation. 

Apart from the above effective, phenomenological 
models, quantum gravitational theories exist that also 
suggest gravitons could have a mass, such as for example 



loop quantum cosmology (LQG) |ll2l . Ilisj . In LQG, the 
dispersion relation of propagating tensor modes acquires 
holonomy cor recti ons during loop quantization, which 
g = A~^/^^~^{p/pc), where 7 is the 



lead to a mass [11 
Barbero-Immirzi parameter, A is related to the area op- 
erator, and p and pc are the total and critical energy den- 
sity respectively. In string theory inspired effective theo- 
ries , the four-dimensional graviton might a lso a cquire a 
mass if there are compact, extra dimensions [ll5| . Unfor- 
tunately, a complete action for such a string theoretical- 
inspired, massive graviton theory that remains covariant 
and free of the van Dam-Veltman-Zakharov (vDVZ) dis- 
continuity (in the limit as nig — *■ 0, the effective theory 
does not reduce to GR) remains elusive. 

Massive graviton theories do not generally satisfy any 
of the criteria described in the Introduction. Some 
of th ese mode ls, including the string-theory inspired 
ones jll5l . Ill6l | , are not weak- field consistent due to the 
vDVZ discontinuity. Moreover, most of these models un- 
avoidingly need to postulate the existence of prior geom- 
etry, making them Lorentz violating and forcing them 
to violate the Strong Equivalence Principle. Finally, the 
main effect of massive graviton theories is to modify GW 
propagation, which is by definition a weak- field effect. 
Needless to say, due to these problems, only limited ef- 
fort has gone into the study of stable solutions [ll7] or 
Gauchy well-posedness, although the latter is unlikely to 
be generically satisfied. 

In spite of such fundamental problems, the GW com- 
munity has been interested in GW propagation tests of 
massive graviton theories because they hold the poten- 
tial to constrain a formally quantum gravitational ef- 
fect. Such tests postulate the phenomenological, special- 
relativistic parameterization [14| : 
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where Eg is the energy of the graviton, Vg its speed and 
the second equality arises in the limit of nig/E ^ 1 . Such 
a modified dispersion relation leads to a modified relation 
between graviton time of emission and arrival [14i] : 
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where fe and /g are the frequencies of the two gravitons 
at emission and D/Dl ^ 1 — z + 0{z'^) in the limit of low 
redshift [l3| . Such a change in the time of arrival leads to 
a modified Fourier transform of the response function [T^ 



(A9) 



where the massive graviton correction is ^'^ 
—PmgU"^ with 
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As before, Ai = rf'/^M is the chirp mass, where M is 
the total mass and 77 is the symmetric mass ratio. 

There are numerous suggestions on how GW detections 
could be used to pl ace bounds on \g [3, [iM [3 [13, fisl . 
[l3,[23,[2il,[iii,[ili]. such tests couM lead to constraints 
as strong as present weak-field bounds, the best of which 
is mg < 3.6 X 10"^^ eV [15^ or equivalently Ag > 3.4 x 
10^^ km [5^ from pulsar timing observations. 

Another feature of a massive graviton that is usually 
neglected is its effect on Newton's second law. A general 
property of massive graviton theories is the introduction 
of a Yukawa-like suppression of the gravitational poten- 
tial fil 



M 



(All) 



Such a modification is sometimes referred to as a fifth 
force, which has led to the Solar System bounds on nig 
cited above p^ . A slightly more general parameteriza- 
tion of the above Yukawa-like correction to the Newto- 
nian potential is 
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where 7mg is an undetermined parameter. Gravitational 
theories with compactified extra dimensions give such 
corrections to the effective 4-dimensional gravitational 
potent ial, w here 7mg depends directly on the compactifi- 
cation jl2l| . Regardless of the Yukawa parameterization, 
GWs will also sense this modified potential as a change in 
the frequency evolution, although the correction is likely 
suppressed by the ratio of the binary separation to the 
graviton Compton wavelength. 

One last important effect in massive graviton theories 
is the ap pear ance of three additional longitudinal modes 
(see eg. |l22l |). Such an effect seems unavoidable if the 
massi ve gr aviton correction to the action is of Pauli-Ficrz 
type [1221 . Longitudinal modes arise due to the non- 
vanishing of the NP scalars 4*2 and ^E's, and can be asso- 
ciated with the presence of spin-1 particles. 

Under the influence of longitudinal GWs, a circular 
arrangement of test particles will be deformed along the 
direction of propagation (see for a graphical repre- 
sentation), which could lead to significant biases in GW 
astronomy. Consider, for example, a situation where a 
GW is observed by a network of detectors, such that 
the inclination angle can be extracted with some accu- 
racy. If GWs were longitudinal instead of transverse, 
the estimation of the inclination angle could have con- 
siderable systematic errors. If there are observed electro- 
magnetic counterparts, for example if binary NS mergers 
cause gamma-ray bursts (GRB), such a bias might be 
noticed with coincident detections. However, if only the 
GRB is observed, and its sky location is used to search 
for a corresponding GW event, the latter could be missed 
if most of the energy was emitted in a longitudinal wave 
not accounted for in the post-analysis. 



3. Chern-Simons Modified Gravity 

CS modified gravity is a particular, non-minimally cou- 
pled, scalar tensor theory, inspired by string theory and 
loop quantum gravity. This extension proposes the in- 
clusion of a parity-violatin g Po ntryagin density to the 
Einstein-Hilbert action (see jl23l | for a complete review) : 

Scs = kJ d'^xy^R+j J d^xy/^^9*RR 
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where *RR — *R" p'^^ R^ a-yS, the dual to the Riemann 
tensor is *R°'p'^^ = fSaxR"^'^^ and k"^ = IGttG. 
The quantities a and (3 are coupling constants, while 
the quantity is a coupling field, and together they 
determine the strength of the CS modification. Note 
that the Pontryagin density is a topological term in four- 
dimensions, and thus, it can usually be integrated by 
parts and removed from the action. In CS modified grav- 
ity, however, this density is coupled to a non-constant 
field, which forces non-trivial correction upon integration 
by parts. A non-dynamical version of this theory exists, 
defined by Eq. (jA13|) with (3 = 0, but this model has been 
shown to be severally overconstraincd [50, 124, 125]. 

The CS modification is motivated by the require- 
ment that any quantum gravitational theory be anomaly- 
free. In perturbative string theory of type I, IIB or 
Heterotic, the absence of such a term leads to Green- 
Schwarz anomalies upon compactification. The require- 
ment of anomaly-cancellat ion thus necessitates the in- 
clusion of such a CS term jl26l |. Even in the non-linear 
regime, Ramond-Ramond scalars induce such a CS term 
for a ll types of string theory due to duality symme- 
tries pa . [1271 

A CS extension of GR can also be motivated by loop 
quantum gravity. Here, such terms arise upon the pro- 
motion of the Barber o-Immirz i parameter to a field when 
coupled to fermions [l28l . Il29l | . Recently, an embedding 
of such a loop quantum gravity effect has been worked 
out in the context of supersymmetric string theory [l3C| . 

Upon variation of the action (jA13|) with respect to all 
degrees of freedom wc find the equations of motion of the 
theory: 
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where T^"°* is the stress-energy for non-gravitational 
matter degrees of freedom, while T^^, is the stress-energy 
tensor of the d field, namely 
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The CS correction is then mostly encoded in the C- 
tensor, which is defined as 

C^"' = [d-y-d) e"'^'^^'V,R''\ + {VjS^) *R^^f"'^^, (A16) 

where e^'"'''' is the Levi-Civita tensor and is the co- 
variant derivative associated with g^^,. 

The main criteria discussed in the Introduction are sat- 
isfied by CS modified gravity. By construction, this GR 
extension is a metric theory that reduces to GR in the 
weak-field limit for sufficiently small gradients of t9. Also, 
as opposed to scalar-tensor theories, CS modified gravity 
can potentially have large, purely strong-field deviations 
due to the presence of the Riemann squared term in the 
action. 

Some of the additional theoretical criteria discussed in 
the Introduction are also satisfied by CS modified grav- 
ity. This GR extension is well-motivated by the leading 
candidates for a quantum theory of gravity. CS modified 
gravity admits the Schwarzschild solution, and a modi- 
fied Kerr solution [13, IsiJ, Il3lj , although the stability of 
these solutions with a dynamical scalar field (with non- 
vanishing background value) have not yet been investi- 
gated. The well-posedness of the theory as a Cauchy 
problem is difficult to determine and it has not yet been 
worked out, although there are arguments that suggests 
the theo ry is well-posed at the linear level (for details on 
this see [l23l |). 

One of the main effects of the CS modification is to in- 
duce parity-violation in gravitational interactions. In the 
non-dynamical framework (when the CS scalar is fixed a 
priori and devoid of dynamics), parity violati on can be 
searched for in the Solar System Il32l.ll33l . Il34| and with 
double binary pulsar observations 135| | through frame- 
dragging modifications. The best current bound on non- 
dynamical CS modified gravity comes from the latter, 
leading to the constraint > 33 mcV. In the dy- 
namical framework (when the CS scalar is consistently 
determined by its coupled evolution equation), the CS 
modification is essentially unconstrained, as the known 
corrections only appear in the strong- field regime [50l.[5l|. 
Moreover, the effects of CS terms can be enh anced in this 
regime in the presence of dense matter |l36l |. 

GW tests have been proposed to const rain CS rn odified 
gravity during binary inspirals [H, HO, Il3ll . Il37j | . Such 
tests can be div ided into those that study GW propaga- 
tion effects [1^ Il37| and th ose that concentrate on GW 
generation effects [SOl Il3l| . GW propagation is mod- 
ified because left- and right-polarized modes obey dif- 
ferent evolution equations, leading to the solution [26| 



hR,L 
given by 



~ "■R,L 



where the CS phase correction is 
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and where Xr^l — ±1 depending on whether the GW is 
right or left-polarized. Note that the correction is purely 
imaginary and fc-dependent, leading to a frequency- 
dependent amplitude correction. The Fourier transform 



of the response fimction of a GW detector is then modi- 
fied as follows ^137j : 
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where i^+,x are the beam-pattern functions of the detec- 
tor, C = cost with b the inclination angle. In Eq. (lAlSp . 
we have effectively neglected terms quadratic in and 
thus, we have ignored the possibility of a CS-indu ced r es- 
onance present in the detector's sensitivity band |137| . 

CS modified gravity also changes the generation of 
GWs, directly affecting their phase and amplitude. Al- 
though these corrections are quadratic in d during the 
early inspiral, they can be large during the late inspiral 
and merger of compact objects. Strong non-linearities 
source large iJ-gradients, which in turn emit energy- 
momentum. Such energy emission modifies the GR flux 
formula, and thus, the orbital frequency evolution, lead- 
ing to phase corrections. Recent studies hav e performed 
simulations of EMRIs in this theory [s^, Il3l| , which con- 
firm that dynamical strong-field effects can lead to sig- 
nificant accumulative dephasing relative to GR. 



4. Einstein- Aether Theory 

Einstein- Aether theory is a Lorentz-violating model 
where the Einstei n-Hilbert action is enlarged into (ssl . 

[T M [1391 [Tiol . [Tin : 
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and ci_2.3,4 are coupling constants. Clearly, the action re- 
mains covariant, but it does not remain Lorentz invariant 
because the vector u^, the "aether," selects a preferred 
frame (eg. that in which is at rest). This vector can 
be forced to have unit norm via an additional Lagrange- 
multiplier term in the action. 

The motivation to construct a theory that explic- 
itly possess preferred-frame physics arises from quan- 
tum gravity. First, the existence of a unique vacuum 
in the quantum theory is believed to select a preferred 
frame 53]. Second, exact Lorentz invariance is known 
to lead to quantum field theoretic divergences for states 
with arbitrarily high energy [138], which are usually elim- 
inated through the introduction of a short-scale cut-off 
that breaks Lorentz invariance. Einstein- Aether theory 
is a model constructed to explicitly break Lorentz invari- 
ance, while preserving the general covariance of GR. 

Einstein- Aether theory satisfies all of the necessary cri- 
teria discussed in the Introduction. As is obvious from 
the form of the action, this model is a metric theory, 
which reduces to GR in the weak field for sufficiently 
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small coupling constants Ci^2,3,4 [l42l Il43l |. In the dy- 
namical strong-field, however, it could disagree with GR, 
although the aether correction to the action is not con- 
structed from higher-order curvature terms. Some of 
the additional theoretical criteria discussed in the In- 
troduction are also satisfied. Although this model is 
phenomenological, it is well-motivated by quantum field 
theory considerations. Stationary and st able BH solu- 
tions have been found both analytically |l44| and nu- 
merically jl45l | . While the well-posedness of the Cauchy 
problem has not been studied separately, the theory has 
been shown to be well -posed numerically in certain dy- 
namical scenarios |l45| . 

Linearized Einstein- Ae ther theory possesses 5 propa- 
gating degrees of freedom [146*1: one scalar (spin-0), three 
vectorial (spin-1) and two tensorial (spin-2) modes. The 
spin-0 mode leads to a breathing mode, while the spin-1 
modes lead to longitudinal polarizations and the spin- 
2 ones to transverse-traceless or quadrupolar ones. The 
speed of propagation of these modes is found to be gener- 
ically different from unity, where in particular for the 
spin-2 modes ^146] 

(A21) 
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In spite of their apparent superluminal nature, the spin- 
2 modes possess positive-definite energies. Although 
Einstein- Aether theory predicts the existence of scalar 
and vectorial modes, these decouple in the weak- field 
limit and all radiation can be forced to be quadrupolar 
by some choice of the coupling constants [111 • 

Com pact objects are also modified in Einstein- Aether 
theory [Til El [US, [HI, [HI [HI . Non-rotating NSs 
have a slightly smaller radii relative to GR, which leads 
to larger redshift factors for electromagnetic radiation 
produced at NS surfaces. The Roche lobes of interacting 
NS is also somewhat enlarged, leading to prompter tidal 
disruptions relative to the GR prediction. BHs possess 
slightly larger inner-most stable circular orbits (ISCO) 
relative to GR, which shifts the GW signal to slightly 
larger f requencie s during the late inspiral, plunge and 
merger jl44l . Il45l | . 

The evolution of binary systems i s also modified in 
Einstein- Aether theory [sl, Il49l Il50l |. Kepler's second 
law is modified in a manner that is dependent on the sen- 
sitivities of the bodies Such a correction should be 
magnified during the late stage of coalescence, precisely 
when the system emits the most amount of energy in GW 
radiation. After the merger, the ringdown phase is also 
modifie d, si nce BHs possess QNR frequencies larger than 
in GR [T49l . [l50t . 



5. MOND and TeVeS 

Modified Newtonian Dynamics (MOND) is a model 
proposed to explain gala xy ro tation curves without in- 
voking cold dark matter [l5ll |. This model was gener- 



alized to a relativistic theory through t he p ostulate of 
a Tensor- Vector-Scalar (TeVeS) theory [nS, EH [HI. 
TeVeS and MOND have managed to successfully explain 
a number of astrophysical phenomena beyond galaxy ro- 
tation curv es, b ut we shall not discuss those here (for a 
review, see [l53j |). 

The main modification introduced in MOND is the cor- 
rection of Newton's second law via 
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where (f> is the Newtonian potential, qq 
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m s is an intrinsic scale in the theory, and fi is 



a function that interpolates between /i(x) ^ x for x 1 
and /i ^ 1 as X ^ oo. One popular choice for this func- 
tion is /i = x{l + x)^^, which leads to the modified equa- 
tion 
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where we have expanded in |VV| 3> ao and we have 
assumed the zeroth order solution a* ~ — V*0. 

How relevant is this modification to binary BH coales- 
cences? The acceleration to compare ao to is |VV| ^ 
M/r^2; where, as before, M is the total mass and ri2 
is the binary separation. The initial separation for a bi- 
nary system with a period such that it is one-year from 
reaching 6M is 
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Supermassive BH binaries (LISA sources) and NS bi- 
naries (ground-detector sources) can be sampled to ini- 
tial separations of ~ 10''Mq and ~ 300Mq before co- 
alescence respectively. For the most massive BH bina- 
ries AI ^ 1O^M0, with the largest initial separations 
ri2 — 1O*M0, the Newtonian acceleration is very low, 
lO^'^ m s~^, but still much larger than ao. Based on 
this, MOND-like modifications to binary dynamics are 
negligibly small. 

A relativistic extension of MOND exists, namely 
TeVeS, which adds a dynamical vector and scalar field 
to the action, making the theory stratified: the gradient 
of the scalar field and the vector field select "strata" or 
preferred spacetime sections. This model is by definition 
bimetric, where the physical metric differs from the aux- 
iliary one, not only by a conformal factor but also by 
terms that depend on the vector field. The interpolating 
function fj, is here governed by the evolution of the scalar 
field, which itself is coupled to the evolution of the met- 
ric tensor. In view of this, one can think of TeVeS as a 
generic theory that interpolates between MOND in the 
slow- veloc ity limit and BD theory in the large-velocity 

limit [m. 

Do these theories satisfy the necessary criteria dis- 
cussed in the Introduction? As we mentioned above. 
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TeVeS is a theory that depends on metric tensors only, 
and as such, it is a metric theory, albeit Lorentz violating. 
In this vein, both MOND and TeVeS violate the Strong 
Equival ence Principle, due to the presence of preferred 
frames [l54| . Weak-field consistency is actually violated 
since TeVeS reduces to MOND in this limit, instead of 
GR or regular Newtonian gravity, though of course if dark 
matter is due to an inaccurate description of gravity on 
galactic scales, then consistency with Newtonian gravity 
here is not a desirable feature. For TeVeS, in the dynam- 
ical strong-field one could expect large deviations from 
GR due to the presence of additional, dynamical scalar 
and vectorial fields. MOND is not well-motivated by fun- 
damental theoretical considerations. TeVeS, being a co- 
variant theory derivable from an action is arguably more 
well motivated, though it is unknown whether i t adm its 
a well-posed initial boundary- value formulation |l55| . 

One consequence for GW detection in TEVES would 
be the modification of the balance law in binaries, since 
additional scalar a nd vectorial degrees of freedom will 
also emit energy [15(^. Another consequence is that 
TeVeS GWs propagate with speed s Vg = c e"*^ < c, 
where (j) is the TeVeS scalar field jlll| . In regions of 
higher curvature (near singularities or close to binary 
mergers) </) is la rger and thus GWs travel more slowly 
than in GR [mj . 



6. DGP Theory 

Several higher-dimensional theories have been pro- 
posed; one ~ Dvali-Gabadadze-Porrati (DGP) pUflSTt- 
has received particular attention in the literature since it 
has the potential to explain the acceleration of the uni- 
verse. Let us then summarize some of the principal fea- 
tures of DGP theory, following the review paper 158,]. 

DGP is a metric theory, based on the string-inspired 
idea of braneworlds: higher-dimensional spaces of which 
our AD universe is only a surface, which for DGP is of co- 
dimension one. The extra dimension in DGP is assumed 
large and flat relative to astrophysical scales. As such, 
the theory introduces infrared modifications to the field 
equations, as opposed to ultraviolet ones, via leakage of 
gravitons into the extra dimension. 

By construction then, DGP gravity will not modify 
GR in the dynamical strong-field regions of relevance to 
compact object GW sources. Nonetheless, it is informa- 
tive to study how such modifications could modify GWs 
emitted by cosmological binaries. Moreover, although 
the theory is somewhat well-motivated from a fundamen- 
tal physics standpoint and although BH and weak-field. 
Solar System solutions exist, it remains unclear whether 
Cauchy well-posedness is satisfied due to the potential 
emergence of g hosts or tachyon-like fields by the extra 
dimension [l59|. 

DGP has not really been studied in the context of GWs 
produced by binary systems. The only studies on GWs 
and DGP refer to either cosmological tensor perturba- 



tions or shock waves [l59|, Il60l Il6l| , neither of which are 
applicable for our purposes. Such a lack of interest in 
GW solutions is surprising, considering the DGP action 
is well-known, in vacuum reducing to 

Sdgp = I d'xV^ (M^i?^^) + Mp^i?^ 

(A25) 

where M and Mp are the five-dimensional and ob- 
served, four-dimensional Planck scale, 7?^^^ is the five- 
dimensional Ricci scalar, (jr(^^^) are the determinants of 
the induced-four and five-dimensional metric tensors and 
Cm is some additional matter lagrangian that couples to 
the four-dimensional sector only. Note that we have rein- 
stated the factors of G here (through the explicit appear- 
ance of the Planck mass), as these get modified in DGP 
theory. With such an action, the DGP field equations 
become, in vacuum 
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where A,B,... — {0, 1, 2, 3, 4} stand for bulk indices and 
. . . — {0, 1, 2, 3} for spacetime indices, with the 
extra dimension, z{x^) the location of the brane and 
^AB^ the four- and five-dimensional Einstein tensor. 

One can now linearize the field equations about 
Minkowski spacetime to find that the spacetime compo- 
nents of the field equations become (in Lorentz gauge): 

Mp{x^)D<^^^h^, ^ Mp{x^)h%^^,-M''D^^^h^,, (A27) 

where □(^'5) are the four- and five-dimensional 
D'Alembertian operators in flat space and we have chosen 
the brane to be located at x^ = 0. If we now re-introduce 
matter fields, which couple to the four-dimensional sec- 
tor of the field equations minimally, we can compute the 
Fourier transform of the above equation, whose solution 
becomes 
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The propagation of GWs is then modified at scales 
larger than ro = Mp/{2M^), which correspond to cosmo- 
logical scales. Alternatively, one can think of DGP mod- 
ifications as introducing a spacetime-dependent Newton 
constant. For example, for scales r < tq the gravitational 
potential of an isolated point source of mass m is 
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where Gtrane = Mp and Gbuik = M 
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One sees that 

such a correction to the binding energy can be mapped 
to the model-independent corrections of Sec. IIII Al via 
p — 1. However, the cross-over scale here is gigantic, as 
since M - lOOMeV then ro ~ IO^^Mq. Of course, if M 
is larger, then the cross-over scale becomes smaller, but 
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then DGP would introduce ultraviolet corrections instead 
of infrared ones. 

An interesting DGP correction to the GW response 
function would occur in the amplitude via the luminos- 
ity distance. Such a correction arises because DGP intro- 
duces modifications to the Friedmann equations and its 
solutions. In particular, the Hubble parameter is modi- 
fied in DGP theory via 
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for a matter dominated and spatially flat cosmology, 
where is the matter density content of the universe. 
Since the luminosity distance depends on an integrated 
measure of the Hubble parameter, this quantity is modi- 
fied in DGP theory. Such a correction could be recovered 
by the ppE parameter a, which models amplitude devi- 
ations, for small redshift sources. 



APPENDIX B: EXOTIC COMPACT OBJECTS 

Many alternative theories of gravity have been dis- 
cussed so far, but little has been said about the generic 
and model-independent possibility that a completely 
gravitationally collapsed object is not represented by the 
Schwarzschild or Kerr line elements. As discussed in the 
Introduction, though many observations of compact ob- 
jects consistent with BHs have been made, there is little 
to no strong evidence that they are indeed BHs as de- 
scribed by GR (ie. that they possess an event horizon). 

Several alternatives exist that replace the 
Schwarzschild or Kerr singularities with other so- 
called exotic alternati ves. One such possibility arises 
from q- balls (see eg. 



162|, [163|), namely a coherent 



scalar "condensate" that can be described classically 
as a non-topological soliton. These objects acquired 
some recent interest since they are unavoidable in 
viable supe rsymmetric extensions of the Standard 
Model [l64| . When the gravitational interaction is 
included, q-stars emerge; furthermore, with gravity, a 
large class of scalar field matter sources admit similar 
stable solitons, an d th e se are of t en c alled boson stars 
(see eg. [165, 16S, EM IHl [l6i [iT^. Boson stars are 
fairly generic objects that could be present not just in 
GR, but also in other a ltern ative theories of gravity, such 
as scalar-tensor ones |l7l| . Some interesting features 
of boson stars are, for example, angul ar mome ntum 
quantization in terms of the scalar charge [I70l . ll72j | and, 
that (even in pure GR) they can be almost as compact 
as BHs, though of course without an event horizon. 



Horizonless, compact objects with large spins are prob- 
lematic from a theoretical standpoint. In fact, boson 
stars, and all horizonless, compact objects that posses 
high spins hav e be en shown to be unstable under small 
perturbations jl73l |. Such instability would lead to grav- 
itational collapse to BHs, or possibly "explosions" . Since 
many astrophysical BH candidates are believed to have 
high spins, such an instability restricts the interest of 
horizonless objects. Nonetheless, the existence of slowly- 
spinning or non-spinning, horizonless, compact objects 
cannot be ruled out by present observations. 

Several works have looked at how merger dynamics 
might be different with boson stars, within the context 
of GR. One study explored the merger of a s mall compact 
object with a supermassive boson star [l74| . and showed 
that stable orbits could exist inside the surface of the 
boson star. Such orbits exist because the effective po- 
tential for these geodesies inside the surface of the boson 
star does not present the usual Schwarzschild-like singu- 
lar behavior, but instead turns over and remains finite, 
allowing for a new minimum. Such orbits lead to extreme 
precession that excite higher frequency harmonics of the 
waveform. 

The merger and ringdown of boson star binaries is 
also drastically different from that of BH binaries. The 
merger of boson stars must be treated numerically; one 
such study explored a scenario where the merger results 
in a single object in a spinning bar configurati on th at ei- 
ther fragments or collapses into a Kerr BH [8^, Il75j . The 
quasinormal ringing of boson stars is also different from 
that of BHs; furthermore the former do not satisfy the 
two-hair (or no-hair) theorem [28^. Instead, boson stars 
satisfy a three-hair theorem, where knowledge of the final 
mass, spin angular momentum and qua drup ole structure 
determines all successive QNR modes |l76j |. A different 
QNR spectrum in GW emission could be used t o sea rch 
for the existence of such boson stars [1^, [2^ [3l|, Il76j | . 

Boson stars, however, are not the only exotic objects in 
the literature. Another interestin g example is a gravita- 
tional vacuum star, or gravastars |177l . Il7a |. In such ob- 
jects, the event horizon is replaced by a phase transition 
that takes the outer Schwarzschild solution to an inner 
De-Sitter spacetime. Once more, the removal of an event 
horizon will lead to modified orbital dynamics, although 
a full analysis has not yet been carried out . Such s tars 
will also lead to a different QNR spectrum |l79l . Il80| | . 

Exotica is not limited to regular, horizonless alter- 
natives; for example, one could have BH solutions in 
GR that have exotic matter "hair" , such as Skyrmion 
BHs jUl. These "dirty" BHs might or might not cause 
modified orbital dynamics, b ut th ey will certainly lead 
to a different QNR spectrum |l82j |. 
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